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an M.Phil. which resulted in writing a FORTRAN program to solve the matrix 
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1.1 Historical review 

Throughout history mankind has tried to control the world in which he lives. From 
the earliest days he realized that his puny strength was no match for the creatures 
around him. Tie could only survive by using his wits and cunning. His major asset 
over all other life forms on earth was his superior intelligence. Stone Age man devised 
tools and weapons from flint, stone and bone and discovered that it was possible to 
train other animals to do his bidding - and so the earliest form of control system was 
conceived. Before long the horse and ox were deployed to undertake a variety of 
tasks, including transport. It took a long time before man learned to replace animals 
with machines. 

Fundamental to any control system is the ability to measure the output of the 
system, and to take corrective action if its value deviates from some desired value. 
This in turn necessitates a sensing device. Man has a number of ‘in-built’ senses 
which from the beginning of time he has used to control his own actions, the actions 
of others, and more recently, the actions of machines. In driving a vehicle for 
example, the most important sense is sight, but hearing and smell can also contribute 
to the driver’s actions. 

The first major step in machine design, which in turn heralded the industrial 
revolution, was the development of the steam engine. A problem that faced engineers 
at the time was how to control the speed of rotation of the engine without human 
intervention. Of the various methods attempted, the most successful was the use of 
a conical pendulum, whose angle of inclination was a function (but not a linear 
function) of the angular velocity of the shaft. This principle was employed by James 
Watt in 1769 in his design of a flyball, or centrifugal speed governor. Thus possibly 
the first system for the automatic control of a machine was born. 

The principle of operation of the Watt governor is shown in Figure 1.1, where 
change in shaft speed will result in a different conical angle of the flyballs. This in 
turn results in linear motion of the sleeve which adjusts the steam mass flow-rate to 
the engine by means of a valve. 

Watt was a practical engineer and did not have much time for theoretical analysis. 
He did, however, observe that under certain conditions the engine appeared to hunt, 
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Steam 




Fig. 1.1 TheWatt centrifugal speed governor. 

where the speed output oscillated about its desired value. The elimination of hunting, 
or as it is more commonly known, instability, is an important feature in the design of 
all control systems. 

In his paper ‘On Governors’, Maxwell (1868) developed the differential equations 
for a governor, linearized about an equilibrium point, and demonstrated that stabil- 
ity of the system depended upon the roots of a characteristic equation having 
negative real parts. The problem of identifying stability criteria for linear systems 
was studied by Hurwitz (1875) and Routh (1905). This was extended to consider the 
stability of nonlinear systems by a Russian mathematician Lyapunov (1893). The 
essential mathematical framework for theoretical analysis was developed by Laplace 
(1749-1827) and Fourier (1758-1830). 

Work on feedback amplifier design at Bell Telephone Laboratories in the 1930s was 
based on the concept of frequency response and backed by the mathematics of complex 
variables. This was discussed by Nyquist (1932) in his paper ‘Regeneration Theory’, 
which described how to determine system stability using frequency domain methods. 
This was extended by Bode (1945) and N ichols during the next 1 5 years to give birth to 
what is still one of the most commonly used control system design methodologies. 

Another important approach to control system design was developed by Evans 
(1948). Based on the work of Maxwell and Routh, Evans, in his Root Locus method, 
designed rules and techniques that allowed the roots of the characteristic equation to 
be displayed in a graphical manner. 
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The advent of digital computers in the 1950s gave rise to the state-space formula- 
tion of differential equations, which, using vector matrix notation, lends itself readily 
to machine computation. The idea of optimum design was first mooted by Wiener 
(1949). The method of dynamic programming was developed by Bellman (1957), at 
about the same time as the maximum principle was discussed by Pontryagin (1962). 
At the first conference of the International Federation of Automatic Control 
(IFAC), Kalman (1960) introduced the dual concept of controllability and observ- 
ability. At the same time Kalman demonstrated that when the system dynamic 
equations are linear and the performance criterion is quadratic (LQ control), then 
the mathematical problem has an explicit solution which provides an optimal control 
law. Also Kalman and Bucy (1961) developed the idea of an optimal filter (Kalman 
filter) which, when combined with an optimal controller, produced linear-quadratic- 
Gaussian (LQG) control. 

The 1980s saw great advances in control theory for the robust design of systems 
with uncertainties in their dynamic characteristics. The work of Athans (1971), 
Safanov (1980), Chiang (1988), Grimble (1988) and others demonstrated how uncer- 
tainty can be modelled and the concept of the FIoo norm and /z-synthesis theory. 

The 1990s has introduced to the control community the concept of intelligent 
control systems. An intelligent machine according to Rzevski (1995) is one that is 
able to achieve a goal or sustained behaviour under conditions of uncertainty. 
Intelligent control theory owes much of its roots to ideas laid down in the field of 
Artificial Intelligence (AI). Artificial Neural Networks (ANNs) are composed of 
many simple computing elements operating in parallel in an attempt to emulate their 
biological counterparts. The theory is based on work undertaken by Hebb (1949), 
Rosenblatt (1961), Kohonen (1987), Widrow-Hoff (1960) and others. The concept of 
fuzzy logic was introduced by Zadeh (1965). This new logic was developed to allow 
computers to model human vagueness. Fuzzy logic controllers, whilst lacking the 
formal rigorous design methodology of other techniques, offer robust control with- 
out the need to model the dynamic behaviour of the system. Workers in the field 
include Mamdani (1976), Sugeno (1985) Sutton (1991) and Tong (1978). 



1.2 Control system fundamentals 
1 .2.1 Concept of a system 

Before discussing the structure of a control system it is necessary to define what is 
meant by a system. Systems mean different things to different people and can include 
purely physical systems such as the machine table of a Computer Numerically 
Controlled (CNC) machine tool or alternatively the procedures necessary for the 
purchase of raw materials together with the control of inventory in a Material 
Requirements Planning (MRP) system. 

However, all systems have certain things in common. They all, for example, 
require inputs and outputs to be specified. In the case of the CNC machine tool 
machine table, the input might be the power to the drive motor, and the outputs 
might be the position, velocity and acceleration of the table. For the MRP system 
inputs would include sales orders and sales forecasts (incorporated in a master 
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V 



Fig. 1.2 The concept of a system. 

production schedule), a bill of materials for component parts and subassemblies, 
inventory records and information relating to capacity requirements planning. Mate- 
rial requirements planning systems generate various output reports that are used in 
planning and managing factory operations. These include order releases, inventory 
status, overdue orders and inventory forecasts. It is necessary to clearly define the 
boundary of a system, together with the inputs and outputs that cross that boundary. 
In general, a system may be defined as a collection of matter, parts, components or 
procedures which are included within some specified boundary as shown in Figure 
1.2. A system may have any number of inputs and outputs. 

In control engineering, the way in which the system outputs respond in changes to 
the system inputs (i.e. the system response) is very important. The control system 
design engineer will attempt to evaluate the system response by determining a 
mathematical model for the system. Knowledge of the system inputs, together with 
the mathematical model, will allow the system outputs to be calculated. 

It is conventional to refer to the system being controlled as the plant, and this, as 
with other elements, is represented by a block diagram. Some inputs, the engineer will 
have direct control over, and can be used to control the plant outputs. These are 
known as control inputs. There are other inputs over which the engineer has no 
control, and these will tend to deflect the plant outputs from their desired values. 
These are called disturbance inputs. 

In the case of the ship shown in Figure 1.3, the rudder and engines are the control 
inputs, whose values can be adjusted to control certain outputs, for example heading 
and forward velocity. The wind, waves and current are disturbance inputs and will 
induce errors in the outputs (called controlled variables) of position, heading and 
forward velocity. In addition, the disturbances will introduce increased ship motion 
(roll, pitch and heave) which again is not desirable. 



Rudder 

Engines 

Wind 

Waves 

Current 




Position 

Forward Velocity 
Heading 
Ship Motion 
(roll, pitch, heave) 



Fig. 1.3 A ship as a dynamic system. 
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Fig. 1.4 Plant inputs and outputs. 

Generally, the relationship between control input, disturbance input, plant and 
controlled variable is shown in Figure 1.4. 



1.2.2 Open-loop systems 

Figure 1.4 represents an open-loop control system and is used for very simple 
applications. The main problem with open-loop control is that the controlled vari- 
able is sensitive to changes in disturbance inputs. So, for example, if a gas fire is 
switched on in a room, and the temperature climbs to 20 °C, it will remain at that 
value unless there is a disturbance. This could be caused by leaving a door to the 
room open, for example. Or alternatively by a change in outside temperature. In 
either case, the internal room temperature will change. For the room temperature to 
remain constant, a mechanism is required to vary the energy output from the gas fire. 



1.2.3 Closed-loop systems 

For a room temperature control system, the first requirement is to detect or sense 
changes in room temperature. The second requirement is to control or vary the energy 
output from the gas fire, if the sensed room temperature is different from the desired 
room temperature. In general, a system that is designed to control the output of a 
plant must contain at least one sensor and controller as shown in Figure 1.5. 



Desired Value 



Forward Path 



Summing 




Fig. 1.5 Closed-loop control system. 
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Figure 1.5 shows the generalized schematic block-diagram for a closed-loop, or 
feedback control system. The controller and plant lie along the forward path, and the 
sensor in the feedback path. The measured value of the plant output is compared at 
the summing point with the desired value. The difference, or error is fed to the 
controller which generates a control signal to drive the plant until its output equals 
the desired value. Such an arrangement is sometimes called an error-actuated system. 



1.3 Examples of control systems 
1.3.1 Room temperature control system 

The physical realization of a system to control room temperature is shown in Figure 
1.6. Flere the output signal from a temperature sensing device such as a thermocouple 
or a resistance thermometer is compared with the desired temperature. Any differ- 
ence or error causes the controller to send a control signal to the gas solenoid valve 
which produces a linear movement of the valve stem, thus adjusting the flow of gas to 
the burner of the gas fire. The desired temperature is usually obtained from manual 
adjustment of a potentiometer. 



Insulation 




Fig. 1.6 Room temperature control system. 
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Temperature 
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(°C) 



Outside 




Fig. 1.7 Block diagram of room temperature control system. 




Introduction to control engineering 7 



A detailed block diagram is shown in Figure 1.7. The physical values of the signals 
around the control loop are shown in brackets. 

Steady conditions will exist when the actual and desired temperatures are the same, 
and the heat input exactly balances the heat loss through the walls of the building. 

The system can operate in two modes: 

(a) Proportional control : Flere the linear movement of the valve stem is proportional to 
the error. This provides a continuous modulation of the heat input to the room 
producing very precise temperature control. This is used for applications where temp- 
erature control, of say better than 1 °C, is required (i.e. hospital operating theatres, 
industrial standards rooms, etc.) where accuracy is more important than cost. 

(b) On-off control'. Also called thermostatic or bang-bang control, the gas valve is 
either fully open or fully closed, i.e. the heater is either on or off. This form of 
control produces an oscillation of about 2 or 3°C of the actual temperature 
about the desired temperature, but is cheap to implement and is used for low-cost 
applications (i.e. domestic heating systems). 



1.3.2 Aircraft elevator control 



In the early days of flight, control surfaces of aircraft were operated by cables 
connected between the control column and the elevators and ailerons. Modern 
high-speed aircraft require power-assisted devices, or servomechanisms to provide 
the large forces necessary to operate the control surfaces. 

Figure 1.8 shows an elevator control system for a high-speed jet. 

Movement of the control column produces a signal from the input angular sensor 
which is compared with the measured elevator angle by the controller which generates 
a control signal proportional to the error. This is fed to an electrohydraulic servovalve 
which generates a spool-valve movement that is proportional to the control signal. 




Fig. 1.8 Elevator control system for a high-speed jet. 
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Fig. 1 .9 Block diagram of elevator control system. 



thus allowing high-pressure fluid to enter the hydraulic cylinder. The pressure differ- 
ence across the piston provides the actuating force to operate the elevator. 

Hydraulic servomechanisms have a good power/weight ratio, and are ideal for 
applications that require large forces to be produced by small and light devices. 

In practice, a ‘feel simulator’ is attached to the control column to allow the pilot to 
sense the magnitude of the aerodynamic forces acting on the control surfaces, thus 
preventing excess loading of the wings and tail-plane. The block diagram for the 
elevator control system is shown in Figure 1.9. 



1.3.3 Computer Numerically Controlled (CNC) machine tool 

Many systems operate under computer control, and Figure 1.10 shows an example of 
a CNC machine tool control system. 

Information relating to the shape of the work-piece and hence the motion of the 
machine table is stored in a computer program. This is relayed in digital format, in a 
sequential form to the controller and is compared with a digital feedback signal from 
the shaft encoder to generate a digital error signal. This is converted to an analogue 



Computer 

Controller 



Machine Table Movement 




Fig. 1.10 Computer numerically controlled machine tool. 
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Fig. 1.11 Block diagram of CNC machine-tool control system. 



control signal which, when amplified, drives a d.c. servomotor. Connected to the 
output shaft of the servomotor (in some cases through a gearbox) is a lead-screw to 
which is attached the machine table, the shaft encoder and a tachogenerator. The 
purpose of this latter device, which produces an analogue signal proportional to 
velocity, is to form an inner, or minor control loop in order to dampen, or stabilize 
the response of the system. 

The block diagram for the CNC machine tool control system is shown in Figure 1.11. 



1.3.4 Ship autopilot control system 

A ship autopilot is designed to maintain a vessel on a set heading while being 
subjected to a series of disturbances such as wind, waves and current as shown in 
Figure 1.3. This method of control is referred to as course-keeping. The autopilot can 
also be used to change course to a new heading, called course-changing. The main 
elements of the autopilot system are shown in Figure 1.12. 

The actual heading is measured by a gyro-compass (or magnetic compass in a 
smaller vessel), and compared with the desired heading, dialled into the autopilot by 
the ship’s master. The autopilot, or controller, computes the demanded rudder angle 
and sends a control signal to the steering gear. The actual rudder angle is monitored 
by a rudder angle sensor and compared with the demanded rudder angle, to form a 
control loop not dissimilar to the elevator control system shown in Figure 1.8. 

The rudder provides a control moment on the hull to drive the actual heading 
towards the desired heading while the wind, waves and current produce moments that 
may help or hinder this action. The block diagram of the system is shown in Figure 1.13. 




Fig. 1.12 Ship autopilot control system. 
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Fig. 1.13 Block diagram of ship autopilot control system. 



1.4 Summary 

In order to design and implement a control system the following essential generic 

elements are required: 

• Knowledge of the desired value : It is necessary to know what it is you are trying to 
control, to what accuracy, and over what range of values. This must be expressed 
in the form of a performance specification. In the physical system this information 
must be converted into a form suitable for the controller to understand (analogue 
or digital signal). 

• Knowledge of the output or actual value : This must be measured by a feedback 
sensor, again in a form suitable for the controller to understand. In addition, the 
sensor must have the necessary resolution and dynamic response so that the 
measured value has the accuracy required from the performance specification. 

• Knowledge of the controlling device : The controller must be able to accept meas- 
urements of desired and actual values and compute a control signal in a suitable 
form to drive an actuating element. Controllers can be a range of devices, including 
mechanical levers, pneumatic elements, analogue or digital circuits or microcomputers. 

• Knowledge of the actuating device : This unit amplifies the control signal and 
provides the ‘effort’ to move the output of the plant towards its desired value. In 
the case of the room temperature control system the actuator is the gas solenoid valve 
and burner, the ‘effort’ being heat input (W). For the ship autopilot system the 
actuator is the steering gear and rudder, the ‘effort’ being turning moment (Nm). 

• Knowledge of the plant : Most control strategies require some knowledge of the 
static and dynamic characteristics of the plant. These can be obtained from 
measurements or from the application of fundamental physical laws, or a com- 
bination of both. 



1.4.1 Control system design 

With all of this knowledge and information available to the control system designer, 
all that is left is to design the system. The first problem to be encountered is that the 
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Fig. 1.14 Steps in the design of a control system. 
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knowledge of the system will be uncertain and incomplete. In particular, the dynamic 
characteristics of the system may change with time (time-variant) and so a fixed 
control strategy will not work. Due to fuel consumption for example, the mass of an 
airliner can be almost half that of its take-off value at the end of a long haul flight. 

Measurements of the controlled variables will be contaminated with electrical 
noise and disturbance effects. Some sensors will provide accurate and reliable data, 
others, because of difficulties in measuring the output variable may produce highly 
random and almost irrelevant information. 

However, there is a standard methodology that can be applied to the design of 
most control systems. The steps in this methodology are shown in Figure 1.14. 

The design of a control system is a mixture of technique and experience. This book 
explains some tried and tested, and some more recent approaches, techniques and 
methods available to the control system designer. Experience, however, only comes 
with time. 




2 



System modelling 



2.1 Mathematical models 

If the dynamic behaviour of a physical system can be represented by an equation, or 
a set of equations, this is referred to as the mathematical model of the system. Such 
models can be constructed from knowledge of the physical characteristics of the 
system, i.e. mass for a mechanical system or resistance for an electrical system. 
Alternatively, a mathematical model may be determined by experimentation, by 
measuring how the system output responds to known inputs. 



2.2 Simple mathematical model of a motor vehicle 



Assume that a mathematical model for a motor vehicle is required, relating the accel- 
erator pedal angle 6 to the forward speed u, a simple mathematical model might be 



II 


(2.1) 


Since u and 6 are functions of time, they are written u(t) and The constant a 

could be calculated if the following vehicle data for engine torque T, wheel traction 
force F, aerodynamic drag D were available 


T = bd(t) 


(2.2) 


II 


(2.3) 


D = d u{t) 


(2.4) 


Now aerodynamic drag D must equal traction force F 





D = F 
d u(t) = cT 
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Forward 

Speed 

u(t) 

(m/s) 




Accelerator angle 6(t) (degrees) 



Fig. 2.1 Vehicle forward speed plotted against accelerator angle. 



from (2.2) 



d u(t) = cbO{t) 



giving 




0(t) 



(2.5) 



Hence the constant for the vehicle is 



a = 




( 2 . 6 ) 



If the constants b, c and cl were not available, then the vehicle model could be 
obtained by measuring the forward speed u(t ) for a number of different accelerator 
angles 0(t) and plotting the results, as shown in Figure 2.1. 

Since Figure 2.1 shows a linear relationship, the value of the vehicle constant a is 
the slope of the line. 



2.3 More complex mathematical models 

Equation (2.1) for the motor vehicle implies that when there is a change in accelerator 
angle, there is an instantaneous change in vehicle forward speed. As all car drivers 
know, it takes time to build up to the new forward speed, so to model the dynamic 
characteristics of the vehicle accurately, this needs to be taken into account. 

Mathematical models that represent the dynamic behaviour of physical systems 
are constructed using differential equations. A more accurate representation of the 
motor vehicle would be 



Here, du/dt is the acceleration of the vehicle. When it travels at constant velocity, this 
term becomes zero. So then 



e-j- +f u = sO(t) 



(2.7) 
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MO = gO(t ) 
u(t)= 

Hence ( gjf ) is again the vehicle constant, or parameter a in equation (2.1) 



(2.8) 



2.3.1 Differential equations with constant coefficients 

In general, consider a system whose output is x(t), whose input is y{t) and contains 
constant coefficients of values a, b, c,..., z. If the dynamics of the system produce a 
first-order differential equation, it would be represented as 

a ^ + bx = cy(i) (2.9) 

dr 

If the system dynamics produced a second-order differential equation, it would be 
represented by 

d 2 x dx 

a -— 5 - + b~ f4- cx = ey(t) ( 2 . 10 ) 

dr- dr 

If the dynamics produce a third-order differential equation, its representation 
would be 



d 3 x , d 2 x dx . .. ..... 

“w +h ifi +c i; +ex = M,) <2in 

Equations (2.9), (2.10) and (2.11) are linear differential equations with constant 
coefficients. Note that the order of the differential equation is the order of the highest 
derivative. Systems described by such equations are called linear systems of the same 
order as the differential equation. For example, equation (2.9) describes a first-order 
linear system, equation ( 2 . 10 ) a second-order linear system and equation ( 2 . 11 ) a 
third-order linear system. 



2.4 Mathematical models of mechanical systems 

Mechanical systems are usually considered to comprise of the linear lumped para- 
meter elements of stiffness, damping and mass. 



2.4.1 Stiffness in mechanical systems 

An elastic element is assumed to produce an extension proportional to the force (or 
torque) applied to it. 

For the translational spring 



Force cx Extension 
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*i(Q *o(0 




P(f) P(t) 



(a) Translational Spring 
Fig. 2.2 Linear elastic elements. 



m e o( f ) 




m 

(b) Rotational Spring 



If A'; (?) > x 0 {t), then 



P(t) = K{x\(t) - x 0 (t)) 



( 2 . 12 ) 



And for the rotational spring 
If 6{{t) > 0 o (t ), then 



Torque oc Twist 



T (0 = - flo(O) ( 2 - 13 ) 

Note that K, the spring stiffness, has units of (N/m) in equation (2.12) and (Nm/rad) 
in equation (2.13). 



2.4.2 Damping in mechanical systems 

A damping element (sometimes called a dashpot) is assumed to produce a velocity 
proportional to the force (or torque) applied to it. 

For the translational damper 

Force oc Velocity 

P(t) = Cv(t) = C^ (2.14) 

And for the rotational damper 

Torque oc Angular velocity 
df) 

T(t) = Coj(t) = (2.15) 



m 

— ► 






► v(t ) 



A 




(a) Translational Damper 



(b) Rotational Damper 



Fig. 2.3 Linear damping elements. 
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Note that C, the damping coefficient, has units of (Ns/m) in equation (2.14) and 
(Nm s/rad) in equation (2.15). 



2.4.3 Mass in mechanical systems 

The force to accelerate a body is the product of its mass and acceleration (Newton’s 
second law). 

For the translational system 



Force <x Acceleration 

P(t) = mail) = m^- = X ° (2.16) 

at at 1 

For the rotational system 

Torque ex Angular acceleration 

m = Ia { t) = I^ = I^ (2.17) 

In equation (2.17) / is the moment of inertia about the rotational axis. 

When analysing mechanical systems, it is usual to identify all external forces by 
the use of a ‘Free-body diagram’, and then apply Newton’s second law of motion in 
the form: 

Y j F = ma for translational systems 
or 

y^M = /a for rotational systems (2.18) 

Example 2. 1 

Find the differential equation relating the displacements x\(t) and x 0 (t) for the 
spring-mass-damper system shown in Figure 2.5. What would be the effect of 
neglecting the mass? 



/ 




Fig. 2.4 Linear mass elements. 
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Fig. 2.5 Spring-mass-damper system. 



K(x \-x Q ) < 



m 



► 




■* 






Fig. 2.6 Free-body diagram for spring-mass-damper system. 



Solution 

Using equations (2.12) and (2.14) the free-body diagram is shown in Figure 2.6. 
From equation (2.18), the equation of motion is 



E 



F x = ma x 



,.r d".\'o d.Y 0 

Kx\ — Kx 0 = m + C 

at - at 



Putting in the form of equation (2.10) 

m + C + Kx o = Kxi(t) (2. 1 9) 

d t- at 

Flence a spring-mass-damper system is a second-order system. 

If the mass is zero then 



K(x ; 



E F * = ° 

Xo) - c ^ = 0 
at 

Kx ; - Kx 0 = C’E 
at 



Hence 

C + Kx o = Kxi(t) 
at 

Thus if the mass is neglected, the system becomes a first-order system. 



( 2 . 20 ) 
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Fig. 2.7 Flywheel in bearings. 



Example 2.2 

A flywheel of moment of inertia I sits in bearings that produce a frictional moment of 
C times the angular velocity u>(t) of the shaft as shown in Figure 2.7. Find the 
differential equation relating the applied torque T(t) and the angular velocity ut(t). 

Solution 

From equation (2.18), the equation of motion is 



= la 

d Ut 

no - c. = i - 

I - 7 -+ Cut = T{t) (2.21) 



Example 2.3 

Figure 2.8 shows a reduction gearbox being driven by a motor that develops a torque 
T m (t). It has a gear reduction ratio of V and the moments of inertia on the motor 
and output shafts are 7 m and I Q , and the respective damping coefficients C m and C 0 . 
Find the differential equation relating the motor torque T m (t) and the output angular 
position 0 o {t). 




Fig. 2.8 Reduction gearbox. 
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d e m d\ 




Fig. 2.9 Free-body diagrams for reduction gearbox. 

Gearbox parameters 

I m = 5 x 10 _6 kgm 2 
I a = 0.01 kgm 2 
C m = 60 x 10~ 6 Nm s/rad 
C 0 = 0. 1 5 Nm s/rad 
w = 50:1 

Solution 

The free-body diagrams for the motor shaft and output shaft are shown in Figure 2.9. 
Equations of Motion are 

(1) Motor shaft 

2 ^ M - Im 

T m (t)-C m ^-aX(t) = I m ^ 
re-arranging the above equation, 

= (2.22) 



E M = / < 



d 2 6» m 
d t 2 



bX(t) - C 0 



ddo 
d t 



= 'o 



d 2 e 0 

d? 2 



(2) Output shaft 
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re-arranging the above equation. 



X{) bV°dt 2+C °dt 



(2.23) 



Equating equations (2.22) and (2.23) 

-(t (t) - I d "' 9m - C — 
a V d/ 2 at 



+ c d0o 
d? 2 d? 



Kinematic relationships 



- = n 0 m (t) = nd 0 (t) 



dO m _ dd 0 

~dt~ ~ n ~dt 



dt 2 



= n 



df- 



Hence 



l T rA T ^ d9 < 

T m (t) w/ m ^ 2 ^ 



= /, 



d z e 0 

dt 2 



giving the differential equation 



d 2 0 o 



de n 



c ° d/ J 



(4 + » 2 / m ) + (Co + n 2 c m ) = »r m (0 



(2.24) 



The terms (7 0 + /r7 m ) and (C 0 + n 2 C m ) are called the equivalent moment of inertia 7 e 
and equivalent damping coefficient C e referred to the output shaft. 

Substituting values gives 

I e = (0.01 + 50 2 x 5 x 10~ 6 ) = 0.0225 kgm 2 
C e = (0.15 + 50 2 x 60 x 10~ 6 ) = 0.3 Nm s/rad 

From equation (2.24) 



0.0225 ^ + 0.3 ^ = 50 T m (t) 
at 2 cu 



(2.25) 



2.5 Mathematical models of electrical systems 

The basic passive elements of electrical systems are resistance, inductance and capa- 
citance as shown in Figure 2.10. 
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MO 



► '(0 

(b) Inductance 

MO | |_C vJS) 

► /(0 

(c) Capacitance 

Fig. 2.10 Passive elements of an electrical system. 

For a resistive element, Ohm’s Law can be written 

(vi (t) - v 2 (t)) = Ri(t) (2.26) 

For an inductive element, the relationship between voltage and current is 

(v,(f) - vzCO) = (2-27) 

For a capacitive element, the electrostatic equation is 

Q(t) = C( V] (t) - vj(0) 

Differentiating both sides with respect to t 

^ = i{ t ) = C^ (v ‘(0 - v 2 (0) (2.28) 

Note that if both sides of equation (2.28) are integrated then 

(vt(0 - V2(0) = idt (2.29) 

Example 2.4 

Find the differential equation relating v\(t) and v 2 (0 for the RC network shown in 
Figure 2.11. 



Solution 

From equations (2.26) and (2.29) 

vi(0 - v 2 (0 = Ri{t) 
nit) =-^J idt 



(2.30) 
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Fig. 2.11 RC network, 
or 

d\h 

C J = /(/) (2.31) 

substituting (2.31) into (2.30) 

n(t) - V2(0 = RC ^ (2.32) 

Equation (2.32) can be expressed as a first-order differential equation 

RC d ^+v 2 = v l (t) (2.33) 

d? 

Example 2.5 

Find the differential equations relating v\(t) and \h(t) for the networks shown in 
Figure 2.12. 



MO 
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(a) 



0(0 +4(0 ^2 







j 


k 1 




0(0 
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(b) 



Fig. 2.12 Electrical networks. 
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Solution for Network (a) Figure 2.12 
From equations (2.26), (2.27) and (2.29) 



d i 

vi (0 - v 2 (t) = Ri(t) + L — 
at 



vi(t) = c I idt 



or 



dv? 



substituting (2.35) into (2.34) 



or 



MO -M0 = /!C^ + id(c^ 



dv ? d“v? 
n(t)-v 2 (t) = RC^ + LC^ 



Equation (2.36) can be expressed as a second-order differential equation 

d 2 v? dv? 

LC d/ + ^ C d7 + V2 = VlW 

Solution for Network (b) Figure 2.12 
System equations 



vi(0 - v 3 (0 = R\(i\(t) + hit)) 



1 f dv? 

v 3 (t) = - / hdt or Ci — = h(t) 



v 3 (0 - v 2 (0 = fl 2 /' 2 (0 



1 /' dv 2 

v 2 (t) = — / / 2 d/ or C 2 = / 2 (0 



From equation (2.40) 

v 3 (0 = Rihit) + v 2 (0 

Substituting for z 2 (t) using equation (2.41) 



(2-34) 



(2.35) 



(2.36) 

(2.37) 

(2.38) 

(2.39) 

(2.40) 

(2.41) 
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Hence from equations (2.42) and (2.39) 

h(t) = C\ — j.R 2 C 2 -^ + TzW 
d 2 v? dr? 

Substituting equations (2.41), (2.42) and (2.43) into equation (2.38) 



f dv 2 ) f d v? 

V,(0 - | R 2 C 2 + V 2 (0 1 = R\ | /? 2 Ci C 2 

which produces the second-order differential equation 



dv 2 dv 2 
Cl d7 +C2 ^7 



RxRiCx C 2 ^ + (Ri Ci +R x C 2 + R 2 C 2 ) ^ + v 2 = vi(t) 



(2.43) 



(2.44) 



2.6 Mathematical models of thermal systems 

It is convenient to consider thermal systems as being analogous to electrical systems 
so that they contain both resistive and capacitive elements. 



2.6.1 Thermal resistance Rj 

Heat flow by conduction is given by Fourier’s Law 

„ KA(e 1 - 0 2 ) 

q t = 2 

The parameters in equation (2.45) are shown in Figure 2.13. They are 
(0 1 — 0 2 ) = Temperature differential (K) 

A = Normal cross sectional area (m 2 ) 

£ = Thickness (m) 

K = Thermal conductivity (W/mK) 

Qj = Heat flow (J/s = W) 




^ L 



(2.45) 



Fig. 2.13 Heat flow through a flat plate. 
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Equation (2.45) can be written in the same form as Ohm’s Law (equation (2.26)) 

mt) - 0 2 (t)) = RjQ T (t) (2.46) 

where Rj is the thermal resistance and is 



Rj — 



i 

KA 



(2.47) 



2.6.2 Thermal capacitance C T 

The heat stored by a body is 

H(t) = mC p 9(t ) (2.48) 

where 

H = Heat (J) 
m = Mass (kg) 

C p = Specific heat at constant 
pressure (J/kgK) 

9 = Temperature rise (K) 

If equation (2.48) is compared with the electrostatic equation 

Q(t) = CMt) (2.49) 

then the thermal capacitance Cj is 

C T = inC p (2.50) 

To obtain the heat flow Qj, equation (2.48) is differentiated with respect to time 



d H 
d t 



= mCp 



d 9 
d t 



or 

df) 

Qt( 0 = Cj 



(2.51) 



(2.52) 



Example 2.6 

Heat flows from a heat source at temperature 9\{t) through a wall having ideal 
thermal resistance Rj to a heat sink at temperature O 2 U) having ideal thermal 
capacitance Cj as shown in Figure 2.14. Find the differential equation relating 
9\{t) and 9 2 {t). 



Qt(0 



( 01 ( 0 - 02 ( 0 ) 

Rj 



Solution 

(1) Wall: From equation (2.46) 



(2.53) 
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Wall 




Fig. 2.14 Heat transfer system. 

(2) Heat sink: From equation (2.52) 

Qi (0 = Ct ? (2 - 54) 

Equating equations (2.53) and (2.54) 

R t t At 

Re-arranging to give the first-order differential equation 

RjCj^+fh = (hit) (2.55) 



2.7 Mathematical models of fluid systems 

Like thermal systems, it is convenient to consider fluid systems as being analogous to 
electrical systems. There is one important difference however, and this is that the 
relationship between pressure and flow-rate for a liquid under turbulent flow condi- 
tions is nonlinear. In order to represent such systems using linear differential equa- 
tions it becomes necessary to linearize the system equations. 



2.7.1 Linearization of nonlinear functions for small 
perturbations 

Consider a nonlinear function Y = f(x) as shown in Figure 2.15. Assume that it is 
necessary to operate in the vicinity of point a on the curve (the operating point) 
whose co-ordinates are X a ,Y a . 

For the small perturbations AX and AT about the operating point a let 



AT = i 
A Y = y 



(2.56) 
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If the slope at the operating point is 



d Y 
dX 



a 



then the approximate linear relationship becomes 



y = 



d7 

dX 



(2.57) 



Example 2.7 

The free-body diagram of a ship is shown in Figure 2.16. It has a mass of 15 x 10 6 kg 
and the propeller produces a thrust of K n times the angular velocity n of the propeller, 
K„ having a value of 1 10 x 10 3 Ns/rad. The hydrodynamic resistance is given by the 
relationship R = C v V 2 , where C v has a value of 10,000 Ns 2 /m 2 . Determine, using 
small perturbation theory, the linear differential equation relating the forward speed 
v(t) and propeller angular velocity n(t) when the forward speed is 7.5 m/s. 



Solution 

Linearize hydrodynamic resistance equation for an operating speed V a of 7.5 m/s. 

R = C v V 2 

§-= 2C ' v 



d R 
dV 

d R 
dV 



= 2 C v V a 

= 2 x 10000 x 7.5 
= C = 150 000 Ns/m 




Fig. 2.15 Linearization of a nonlinear function. 
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Fig. 2.16 Free-body diagram of ship. 



Hence the linear relationship is 

R = 

Using Newton’s second law of motion 

T — R = 

K„n — Cv = m^- 
at 

dv 

m — + Cv = K„n 
d t 



Cv 

ma x 

dv 
in — 



(2.58) 



(2.59) 



Substituting values gives 

(15 x 10 6 )^ + (150 x 10 3 )v = (110 x 10>(?) (2.60) 

d t 

Example 2.8 

In Figure 2.17 the tank of water has a cross-sectional area A, and under steady 
conditions both the outflow and inflow is V a and the head is H a . 



(a) Under these conditions find an expression for the linearized valve resistance R[ 
given that flow through the valve is 



V = A x CA\flgH . 



where 

V = volumetric flow-rate (m 3 /s) 

A v = valve flow area (m 2 ) 

Cd = coefficient of discharge 
g = acceleration due to gravity (m/s 2 ) 
H = head across the valve (m) 
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(b) If the steady value of the head H a is 1.5 m, what is the valve resistance R\ when 

A v = 15 x 1(T 3 m 2 
g = 9.81 m/s 2 
C d = 0.6 

(c) If the inflow now increases an amount vi producing an increase in head h and an 
increase in outflow V 2 , find the differential equation relating v\ and V 2 when the 
tank cross-sectional area A is 0.75m 2 . 

Solution 

(a) Flow through the valve is given by 



V=A v C d ^2gH 



now 




The linearized relationship is 



h = R[V2 



hence 




(2.61) 



h 



//|\ x 



H. 



A 



s 



Ft, 




Fig. 2.17 Tank and valve system. 
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(b) Inserting values gives 



1 



Rf 15 x 10- 3 x 0.6 V 9.81 
R[ = 61.45 s/m 2 



2 x 1.5 



(c) Tank (Continuity Equation) 



dh 

Inflow — Outflow = A — 
d t 

{V a + Vi) — (V a + V 2 ) = A — 

d t 

dh 

vi - v 2 = A — 
dt 



Valve (Linearized Equation) 



and 



h = R(V 2 



dh n dv ’2 
~r = R[ —7— 
dt dt 



Substituting equation (2.64) into equation (2.63) 

, „ dV2 
vi - v 2 = ARf ~~~ 
dt 



giving 



Inserting values gives 



dv-> 

AR { -^+v 2 = v\(t) 



46.09 + V 2 = vi (?) 

dt 



(2.62) 



(2.63) 



(2.64) 



(2.65) 



(2.66) 



2.8 Further problems 

Example 2.9 

A solenoid valve is shown in Figure 2.18. The coil has an electrical resistance of 4fi, 
an inductance of 0.6 H and produces an electromagnetic force F c (t) of K c times the 
current i(t). The valve has a mass of 0.125 kg and the linear bearings produce 
a resistive force of C times the velocity u(t). The values of K c and C are 0.4 N/A 
and 0.25 Ns/m respectively. Develop the differential equations relating the voltage 
v(t) and current i(t) for the electrical circuit, and also for the current i(t) and velocity 
u(t) for the mechanical elements. Hence deduce the overall differential equation 
relating the input voltage v(t) to the output velocity u(t). 
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FM 



Fig. 2.18 Solenoid valve. 



Solution 



L t, +Ri = ,{ 0 

d u 

m — + Cu = K c i(t) 
at 



d 1 u d u 

0.075 — + 0.65 — + u = 0.4v(0 
dr 2 d t 



Example 2.10 

The laser-guided missile shown in Figure 2.19 has a pitch moment of inertia of 
90kgm 2 . The control fins produce a moment about the pitch mass centre of 
360 Nm per radian of fin angle /3ft). The fin positional control system is described 
by the differential equation 

°-2^y + /3(t) = u(t) 

where u(t) is the control signal. Determine the differential equation relating the 
control signal u(t) and the pitch angle 9{t). 



Solution 



d 3 <9 d 2 9 

di^ +5 dfi 



20u(t) 




Fig. 2.19 Laser-guided missile. 
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Fig. 2.20 Torsional spring-mass-damper system. 



Example 2.11 

A torsional spring of stiffness K, a mass of moment of inertia I and a fluid damper 
with damping coefficient C are connected together as shown in Figure 2.20. If the 
angular displacement of the free end of the spring is 9 x {t) and the angular displace- 
ment of the mass and damper is 9 a (t), find the differential equation relating 9\{t) and 
9 0 (t) given that 

I = 2.5 kg m 2 
C = 12.5 Nm s/rad 
K = 250 Nm/rad 



Solution 

a2q aq 

2.5—^+ 12.5— -+25O0 o = 25O0i(O 
at 1 at 

Example 2.12 

A field controlled d.c. motor develops a torque T m {t) proportional to the field current 
if(t). The rotating parts have a moment of inertia / of 1.5 kg m 2 and a viscous 
damping coefficient C of 0.5 Nm s/rad. 

When a current of 1.0 A is passed through the field coil, the shaft finally settles 
down to a steady speed u Q (t) of 5 rad/s. 

(a) Determine the differential equations relating i[(t) and uj 0 (t). 

(b) What is the value of the coil constant K c , and hence what is the torque developed 
by the motor when a current of 0.5 A flows through the field coil? 

Solution 

(a) /^ + Cw 0 = K c i f (t) 

at 

(b) K c = 2.5 Nm/A. T m = 1 .25 Nm 
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MO 



C 



MO 



nr 



Fig. 2.21 Passive RC network. 



Oven 



Q(0 




* 



K 



► 



Burner 



0,(0 



Fig. 2.22 Drying oven. 

Example 2.13 

Figure 2.21 shows a passive electrical network. Determine the differential equation 
relating vi(?) and V 2 (?)• 



Example 2.14 

A drying oven which is constructed of firebrick walls is heated by an electrically 
operated gas burner as shown in Figure 2.22. The system variables and constants are 

v\{t) = burner operating voltage (V) 

Qi(t) = heat input to oven (W) 

9 0 (t) = internal oven temperature (K) 

9 s (t) = temperature of surroundings (K) 

K = burner constant = 2000W/V 
Rj = thermal resistance of walls = 0.5 x 10 _3 min K/J 
Cj = oven thermal capacitance = 1 x 10 4 J/K 

Find the differential equation relating v,(t), 0 o (t ) and 9 s (t). 



Solution 




Solution 



3 



Time domain analysis 



3.1 Introduction 



The manner in which a dynamic system responds to an input, expressed as a function 
of time, is called the time response. The theoretical evaluation of this response is said 
to be undertaken in the time domain, and is referred to as time domain analysis. It is 
possible to compute the time response of a system if the following is known: 

• the nature of the input(s), expressed as a function of time 

• the mathematical model of the system. 

The time response of any system has two components: 

(a) Transient response : This component of the response will (for a stable system) 
decay, usually exponentially, to zero as time increases. It is a function only of the 
system dynamics, and is independent of the input quantity. 

(b) Steady-state response : This is the response of the system after the transient 
component has decayed and is a function of both the system dynamics and the 
input quantity. 




Fig. 3.1 Transient and steady-state periods of time response. 
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The total response of the system is always the sum of the transient and steady-state 
components. Figure 3.1 shows the transient and steady-state periods of time 
response. Differences between the input function x,(t) (in this case a ramp function) 
and system response x Q (t) are called transient errors during the transient period, and 
steady-state errors during the steady-state period. One of the major objectives of 
control system design is to minimize these errors. 



3.2 Laplace transforms 

In order to compute the time response of a dynamic system, it is necessary to solve 
the differential equations (system mathematical model) for given inputs. There are 
a number of analytical and numerical techniques available to do this, but the one 
favoured by control engineers is the use of the Laplace transform. 

This technique transforms the problem from the time (or t) domain to the Laplace 
(or s ) domain. The advantage in doing this is that complex time domain differential 
equations become relatively simple s domain algebraic equations. When a suitable 
solution is arrived at, it is inverse transformed back to the time domain. The process 
is shown in Figure 3.2. 

The Laplace transform of a function of time/(t) is given by the integral 




(3.1) 



where 5 is a complex variable a ± jco and is called the Laplace operator. 



s Domain F(s) 



Laplace 

Transform 



> Algebraic 
equations 



5W)] = F(s) 



Inverse 

Laplace 

Transform 



ar\F( S )]=f(t) 



Time Domain f(t) 

Differential 

equations 



Fig. 3.2 The Laplace transform process. 
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3.2.1 Laplace transforms of common functions 

Example 3.1 

f{t ) = 1 (called a unit step function). 



Solution 

From equation (3.1) 



Example 3.2 



,r|/in| = Fivi= / ie-"d i 






Jo 



--( 0 - 1 ) 

5 



fit) = e- 



W(0] = *W = f 

J o 



e~ a, e-' r d r 



-(s+a)t df 



(3.2) 



1 



5 + 0 

1 

5 + 0 



(e+" +a)r ) 

( 0 - 1 ) 



J o 



1 



5 + 0 



(3.3) 



Table 3.1 gives further Laplace transforms of common functions (called Laplace 
transform pairs). 



3.2.2 Properties of the Laplace transform 

(a) Derivatives: The Laplace transform of a time derivative is 

= s n F(s) -/(OK 1 - 1 -/'(0)5"- 2 - • • • (3.4) 

where /(0),/'(0) are the initial conditions, or the values of fit), d/d t fit) etc. at t = 0 

(b) Linearity 



n.m ±fiit)\ = Fiis) ± F 2 is ) 



(3.5) 
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Table 3.1 Common Laplace transform pairs 



Time function f(t) 


Laplace transform £T[f{t)\ — F(s) 


i 


unit impulse S(t) 


i 


2 


unit step 1 


l/s 


3 


unit ramp t 


1// 


4 


f 


n\ 

S n + 1 


5 


e~ al 


1 

(s + a) 


6 


1 - e~“' 


a 

s(s + a) 


7 


sine ot 


CJ 

s 2 +u 2 


8 


coso ;/ 


s 

S 2 +CC 2 


9 


e _< " sin ut 


u 

( s + a f + J 2 


10 


e~ a '(cos ut — — sin uit) 

U J 


s 

(s + a) 2 + U) 2 



(c) Constant multiplication 



JF[af{t)\ = aF(s) 



(d) Real shift theorem 



J F[f(t - D] = e~ Ts F(s) for T > 0 



(e) Convolution integral 



~ T)d t = F { (s)F 2 (s) 



(f) Initial value theorem 

/(0) = I ini [/(/)] = lim [sF(s)] 

t—>0 s — >oo 

(g) Final value theorem 

/( oo) = lim [.AO] = lim [sF(s)] 

t—> oo s— >0 



3.2.3 Inverse transformation 

The inverse transform of a function of 5 is given by the integral 

i ra+ )u) 

f(t) = J?- l [F(s)\ = — F(s)c sl ds 
27TJ Ja- p, 



(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 




Time domain analysis 39 



In practice, inverse transformation is most easily achieved by using partial fractions 
to break down solutions into standard components, and then use tables of Laplace 
transform pairs, as given in Table 3.1. 



3.2.4 Common partial fraction expansions 

(i) Factored roots 

K A B 

s(s + a) s (s + a) 

(ii) Repeated roots 

K A B C 

s 2 {s + a) s s 2 + (s + a) 

(iii) Second-order real roots ( b 2 > 4 ac) 



(3.12) 



(3.13) 



K 

s(as 2 + bs + c) 



K 

s(s + d)(s + e ) 



A B C 

s + (s + d ) + (s + e) 



(iv) Second-order complex roots ( b 2 < 4 ac) 

K _ A Bs + C 
s(as 2 + bs + c) s as 2 + bs + c 

Completing the square gives 

A Bs+C 
s (s + a) 2 + co 2 

Note: In (iii) and (iv) the coefficient a is usually factored to a unity value. 



(3.14) 



3.3 Transfer functions 

A transfer function is the Laplace transform of a differential equation with zero 
initial conditions. It is a very easy way to transform from the time to the s domain, 
and a powerful tool for the control engineer. 

Example 3.3 

Find the Laplace transform of the following differential equation given: 

(a) initial conditions x a = 4, dx 0 /dt = 3 

(b) zero initial conditions 



d 2 



*0 , dxp 

d t 2 d t 



2x 0 = 5 
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*i(s) 



XJ.s) 



► 



G(s) 



► 



Fig. 3.3 The transfer function approach. 

Solution 

(a) Including initial conditions : Take Laplace transforms (equation (3.4), Table 3.1). 



Example 3.3(b) is easily solved using transfer functions. Figure 3.3 shows the general 
approach. In Figure 3.3 

• T,(,v) is the Laplace transform of the input function. 

• A 0 (s) is the Laplace transform of the output function, or system response. 

• G(s) is the transfer function, i.e. the Laplace transform of the differential equation 
for zero initial conditions. 

The solution is therefore given by 




(■ s 2 + 3s + 2)A 0 (,s) 



5 + As 1 + 1 5s 




(3.15) 



(b) Zero initial conditions 
At t = 0, x 0 = 0, dx a /dt = 0. 



Take Laplace transforms 



s L X 0 (s) + 3LL 0 (y> + 2 




(3.16) 



X 0 is) = G{s)X l (s) 



(3.17) 



Thus, for a general second-order transfer function 




(as 2 + bs + c)Ao(^) = KX\(s) 



Hence 
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m 


K 


Us) 




as+bs+c 





Fig. 3.4 General second-order transfer function. 



X-(s)=bts 


1 


Us) 




S 2 +3s+2 





Fig. 3.5 Example 3.3(b) expressed as a transfer function. 



Comparing equations (3.17) and (3.18), the transfer function G(s) is 



G(s) 



K 

as 2 + bs + c 



(3.19) 



which, using the form shown in Figure 3.3, can be expressed as shown in Figure 3.4. 

Returning to Example 3.3(b), the solution, using the transfer function approach is 
shown in Figure 3.5. From Figure 3.5 



X 0 (s) 



5 

5(5 2 + 3^+2) 



(3.20) 



which is the same as equation (3.16). 



3.4 Common time domain input functions 
3.4.1 The impulse function 

An impulse is a pulse with a width At — > 0 as shown in Figure 3.6. The strength of an 
impulse is its area A, where 

A = height h x At. (3.21) 

The Laplace transform of an impulse function is equal to the area of the function. 
The impulse function whose area is unity is called a unit impulse 6(t). 



3.4.2 The step function 

A step function is described as x\(t) = B\ Xi(s) = B/s for t > 0 (Figure 3.7). For a unit 
step function x\(t) = 1; X,(.s) = 1 Is. This is sometimes referred to as a ‘constant 
position’ input. 
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Impulse 




Fig. 3.6 The impulse function. 



*i(?) 



B 



t 



Fig. 3.7 The step function. 



3.4.3 The ramp function 

A ramp function is described as xi(t) = Qt; X[{s ) = Q/s 2 for t > 0 (Figure 3.8). For a 
unit ramp function x\(t) = t; X\(s) = l/s 2 . This is sometimes referred to as a ‘constant 
velocity’ input. 



3.4.4 The parabolic function 

A parabolic function is described as x;(f) = Kt 2 -, T)(,s) = 2 K/s 2 for t > 0 (Figure 3.9). 
For a unit parabolic function Xi(f) = t 2 ; T)(,s) = 2/s 2 . This is sometimes referred to as 
a ‘constant acceleration’ input. 
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t 



Fig. 3.8 The ramp function. 




Fig. 3.9 The parabolic function. 



3.5 Time domain response of first-order systems 
3.5.1 Standard form 

Consider a first-order differential equation 

a + bx Q = cx{{t) (3.22) 

dr 

Take Laplace transforms, zero initial conditions 



asX 0 (s ) + bX Q (s) = cXi(s) 
(as + b)X Q (s) = cX\(s) 
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The transfer function is 



G(s) 



Xo 



(s) 



c 

as + b 



To obtain the standard form, divide by b 



G(s) 



C 

b 



1+5* 



which is written 



G(s) 



K 

1 + Ts 



(3.23) 



Equation (3.23) is the standard form of transfer function for a first-order system, 
where K = steady-state gain constant and T = time constant (seconds). 



3.5.2 Impulse response of first-order systems 

Example 3.4 (See also Appendix 1, examp34.m ) 

Find an expression for the response of a first-order system to an impulse function of 
area A. 



Solution 

From Figure 3.10 



or 



, AK AK/T 
Xo( ' V) " 1 + Ts ~ {s+l/T) 



X 0 {s) 



T \(s + a)) 



(3.24) 



(3.25) 



Equation (3.25) is in the form given in Laplace transform pair 5, Table 3.1, so the 
inverse transform becomes 




—at 




(3.26) 



The impulse response function, equation (3.26) is shown in Figure 3.11. 



X|(s)=/4 



A 


K 


X 0 (s) 




1 + Ts 





Fig. 3.10 Impulse response of a first-order system. 
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*o(f) 




Fig. 3.11 Response of a first-order system to an impulse function of area A. 



3.5.3 Step response of first-order systems 

Example 3.5 (See also Appendix 1, examp35.ni) 

Find an expression for the response of a first-order system to a step function of 
height B. 

Solution 

From Figure 3.12 



Equation (3.27) is in the form given in Laplace transform pair 6 Table 3.1, so the 
inverse transform becomes 



When time t is expressed as a ratio of time constant T, then Table 3.2 and Figure 3.13 
can be constructed. 

Table 3.2 Unit step response of a first-order system 




(3.27) 




(3.28) 



If B = 1 (unit step) and K = 1 (unity gain) then 




(3.29) 



t/T 0 0.25 0.5 0.75 1 

x a (t) 0 0.221 0.393 0.527 0.632 



1.5 2 2.5 3 4 

0.770 0.865 0.920 0.950 0.980 
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X,(s)=B/s 



* 



K 

1 + Ts 



* * 0 (s) 



Fig. 3.12 Step response of a first-order system. 

3.5.4 Experimental determination of system time constant 
using step response 

Method one\ The system time constant is the time the system takes to reach 63.2% of 
its final value (see Table 3.2). 

Method two : The system time constant is the intersection of the slope at t = 0 with 
the final value line (see Figure 3.13) since 

x 0 (t ) = 1 - e~ r/T 

< 3 -»> 

d ^° U = l att = 0 (3.31) 

This also applies to any other tangent, see Figure 3.13. 




Fig. 3.13 Unit step response of a first-order system. 
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3.5.5 Ramp response of first-order systems 

Example 3.6 

Find an expression for the response of a first-order system to a ramp function of 
slope Q. 



Solution 

From Figure 3.14 



X 0 (s) = 



QK 

5 2 (1 + Ts) 



QK/T 

s 2 (s+l/T) 



A | B ( C 
7 + 7 + Q+ \/T) 



(3.32) 



(See partial fraction expansion equation (3.13)). Multiplying both sides by 
s 2 (s + 1 /T), we get 



QK 



= 4di+- 






Cs 2 



i.e. 



QK A 2 A „ B „ 2 
= As 2 + -s + Bs + -+ Cs~ 



Equating coefficients on both sides of equation (3.33) 



From (3.34) 

From (3.36) 

Substituting into (3.35) 
FLence from (3.34) 



(r): 

( 7 ): 

(/): 



0 = A + C 
o = y + b 

QK B 

""t” _ T 



C = —A 



B = QK 



yf = -QKT 



C= QKT 



(3.33) 

(3.34) 

(3.35) 

(3.36) 



7((s)=0/s 2 


K 






1 + Ts 





Fig. 3.14 Ramp response of a first-order system (see also Figure A1.1). 
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Number of Time Constants 



Fig. 3.15 Unit ramp response of a first-order system. 



Inserting values of A, B and C into (3.32) 



Xo(s) 



QKT | QK | QKT 

— V + + (s+l/f) 



Inverse transform, and factor out KQ 



(3.37) 



x 0 (t) = KQ(t-T+ Te^ 7 ) (3.38) 

If Q = 1 (unit ramp) and K = 1 (unity gain) then 

x a (t) = t- T+ 77T' /r (3.39) 

The first term in equation (3.39) represents the input quantity, the second is the 
steady-state error and the third is the transient component. When time t is expressed 
as a ratio of time constant T, then Table 3.3 and Figure 3.15 can be constructed. In 
Figure 3.15 the distance along the time axis between the input and output, in the 
steady-state, is the time constant. 



Table 3.3 Unit ramp response of a first-order system 



t/T 0 1 2 3 4 5 6 7 

Xi (t)/T 0 1 2 3 4 5 6 7 

x a (t)/T 0 0.368 1.135 2.05 3.018 4.007 5 6 
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3.6 Time domain response of second-order systems 
3.6.1 Standard form 

Consider a second-order differential equation 

d 2 .\'r, d.\' n 

a—rr + b—, — I- cx Q = ex\(t) (3.40) 

d?- d? 

Take Laplace transforms, zero initial conditions 

as 2 X 0 (s) + bsX 0 (s) + cX 0 (s) = eX\ (s) 

(as 2 + bs + c)X 0 (s ) = eXi(s) (3.41) 



The transfer function is 



G(s) = ^(s) = 



Xi ‘ as 1 + bs + c 
To obtain the standard form, divide by c 



G(s) = 



-s 2 + -s + 1 



which is written as 



G(s) 



K 



Xs 1 + 1 



This can also be normalized to make the s 2 coefficient unity, i.e. 



G(s) 



Klo 



2 

n 



s 2 + 2 (ui n s + cJ 2 



(3.42) 



(3.43) 



Equations (3.42) and (3.43) are the standard forms of transfer functions for a second- 
order system, where K = steady-state gain constant, w n = undamped natural 
frequency (rad/s) and £ = damping ratio. The meaning of the parameters w n and ( 
are explained in sections 3.6.4 and 3.6.3. 



3.6.2 Roots of the characteristic equation and their 

relationship to damping in second-order systems 

As discussed in Section 3.1, the transient response of a system is independent of the 
input. Thus for transient response analysis, the system input can be considered to be 
zero, and equation (3.41) can be written as 

(as 2 + bs + c)X 0 (s ) = 0 

If A 0 (.y) ^ 0, then 



as 1 + bs + c = 0 



(3.44) 
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Table 3.4 Transient behaviour of a second-order system 



Discriminant 


Roots 


Transient response type 


b 1 > 4 ac 


.si and S2 real 


Overdamped 




and unequal 


Transient 




(-ve) 


Response 


b 2 = 4 ac 


and S2 real 


Critically 




and equal 


Damped Transient 




(-ve) 


Response 


b 2 < 4 ac 


and S2 complex 


Underdamped 




conjugate of the 


Transient 




form: si, S2 = — cr±j uj 


Response 



This polynomial in s is called the Characteristic Equation and its roots will determine 
the system transient response. Their values are 



s\, s 2 = 



—b ± Vb 1 — 4 ac 
2 a 



(3.45) 



The term ( b 2 — 4 ac), called the discriminant, may be positive, zero or negative which 
will make the roots real and unequal, real and equal or complex. This gives rise to the 
three different types of transient response described in Table 3.4. 

The transient response of a second-order system is given by the general solution 

x Q (t) = Ae s '‘ + Be S2 ‘ (3.46) 

This gives a step response function of the form shown in Figure 3.16. 




Fig. 3.16 Effect that roots of the characteristic equation have on the damping of a second-order system. 



Time domain analysis 51 



3.6.3 Critical damping and damping ratio 



Critical damping 

When the damping coefficient C of a second-order system has its critical value C c , the 
system, when disturbed, will reach its steady-state value in the minimum time without 
overshoot. As indicated in Table 3.4, this is when the roots of the Characteristic 
Equation have equal negative real roots. 

Damping ratio £ 

The ratio of the damping coefficient C in a second-order system compared with the 
value of the damping coefficient C c required for critical damping is called the 
Damping Ratio ( (Zeta). Hence 

C = ^ (3.47) 

Cc 

Thus 

£ = 0 No damping 
( < 1 Underdamping 
( = 1 Critical damping 
( > 1 Overdamping 



Example 3.7 

Find the value of the critical damping coefficient C c in terms of K and m for the 
spring-mass-damper system shown in Figure 3.17. 




Lumped Parameter Diagram 
(a) 



Kx 0 Cx 0 



i 






m 




r 



F(t) 



X 0 (!) 
*o(t) 



Free-Body Diagram 
(b) 



Fig. 3.17 Spring-mass-damper system. 
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Solution 

From Newton’s second law 

Fx = mx Q 

From the free-body diagram 

F(t) — Kx 0 (t) — Cx Q (t) = mx 0 (t ) 
Taking Laplace transforms, zero initial conditions 

F(s ) — KX 0 (s ) — CsX 0 (s) = ms 2 X 0 (s) 



Characteristic Equation is 



and the roots are 



(ms 2 + Cs + K)X 0 (s) = F(s) 



ms- + Cs + K = 0 



2 C K n 

i.e. s z H 1 — = 0 

m m 



4- 



For critical damping, the discriminant is zero, hence the roots become 



Also, for critical damping 



5l = S2 = - : 



C: 4 K 



giving 



C c = 2 VKm 



3.6.4 Generalized second-order system response to a unit step 
input 

Consider a second-order system whose steady-state gain is K, undamped natural 
frequency is u> n and whose damping ratio is (, where £ < 1. For a unit step input, the 
block diagram is as shown in Figure 3.18. From Figure 3.18 

Kur 

Xo(s) = , , , ° 

s(s- + 2C,ijJ a S + LU^) 



(3.52) 
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x;(s)=i/s 




Fig. 3.18 Step response of a generalized second-order system for £ < 
Expanding equation (3.52) using partial fractions 



s (s 2 + 2(u> n s + w 2 ) 
Equating (3.52) and (3.53) and multiply by s(s 2 + 2(to n s + w 2 ) 



Equating coefficients 



AC' 2 — A (,v 2 T 2(uj n s ~h u.’ 2 j T Bs~ T Cs 

(s 2 ) : 0 = A + B 

C? 1 ): 0 = 2<w„,4 + C 

(/) : AC 2 = ur n A 



giving 



A = K, B = —K and C = -2(u; n K 
Substituting back into equation (3.53) 

r„(j) = Jc ri 1 s + 2< “" 

Completing the square 

X 0 (s) = K 



1 



s I .s 2 + 2 (w n s + w 2 
5 + 2C > u n 



S I (.S + (Mi) + — C w n 



= A" 



s + 2£w n 



C? + £<U n )“ + (oJn \/ 1 ~ C 2 ) 



(3.53) 



(3.54) 



The terms in the brackets { } can be written in the standard forms 10 and 9 in 
Table 3.1. 



Term (1) = 



—s 



(S + (Mi) + (w a \/ 1 — C 2 ) 



Term (2) = — < 



2(Mi 






"" VT^C-j |(,2 + CWn) 2 + ^ nV ^2y 
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Inverse transform 



K 0 (t) = K 1 - | e Cw "' ^cos (u n y/l - C 2 ) t - = sin (w n \/l - C 2 ) ^ | 

- { TTrbil { e_< " ( sin (“- ' /T3? ) ') } 



Equation (3.55) can be simplified to give 



1 -e <w " ' | cos (uj n \J 1 — f 2 ) t + |-^==j shinny 7 ! - C 2 )?| (3-56) 



Xo(t) = K 1 - 



When C = 0 

x 0 (t) = 7f[l - e°{cos<u n ? + 0}] 

= K[\ -cosw„t] (3.57) 

From equation (3.57) it can be seen that when there is no damping, a step input will 
cause the system to oscillate continuously at w n (rad/s). 

Damped natural frequency to d 

From equation (3.56), when 0 < C> 1, the frequency of transient oscillation is 
given by 

w d = cWl -C 2 (3-58) 

where lo& is called the damped natural frequency. Flence equation (3.56) can be 
written as 

x 0 (t) = K 1 — e~ CWn '|cos w ( i/ + - ^===^ sinwdt| (3.59) 



where 



= K 1 7 sinpudf + <b) 



tan d> = ■ 



When £ = 1, the unit step response is 

x 0 (t) = K[ 1 - e” w "'(l + u D t)\ 

and when £ > 1, the unit step response from equation (3.46) is given by 



x 0 (t) = K 1 - < - + 






+ |i £_) e (-t-V5=i)«. 
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“ n t{ rad) 



Fig. 3.19 Unit step response of a second-order system. 

The generalized second-order system response to a unit step input is shown in Figure 
3.19 for the condition K = 1 (see also Appendix 1, sec_ord.m). 



3.7 Step response analysis and performance 
specification 

3.7.1 Step response analysis 

It is possible to identify the mathematical model of an underdamped second-order 
system from its step response function. 

Consider a unity-gain (K = 1) second-order underdamped system responding to 
an input of the form 

x\(t) = B (3.64) 

The resulting output x Q (t ) would be as shown in Figure 3.20. There are two methods 
for calculating the damping ratio. 

Method (a): Percentage Overshoot of first peak 

%Overshoot = ^ x 100 (3.65) 

B 

Now 

a x = fi e -^" (T / 2) 
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x 0 (t) 



B 



B 




Fig. 3.20 Step response analysis. 



Thus, 



5 e -(“ n(r/2) 

%Overshoot = x 100 (3.66) 



Since the frequency of transient oscillation is u>d, then, 

2n 

UJd 

2n 

UnV 1 - C 2 

Substituting (3.67) into (3.66) 

% Overshoot = x 10 o 

% Overshoot = 



(3.67) 



(3.68) 



Method (b): Logarithmic decrement. Consider the ratio of successive peaks ci\ and a 2 

a\ = Be-^12) (3. 69 ) 

a 2 = Se' f “” (3r/2) (3.70) 

Hence 

_ _ e p {-C“n(-r/2)+Cw„(3r/2)} 

a2 e -C^„(3r/2) 

— gC“nT _ ^n/y/l-C 2 



( 3 . 71 ) 
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Equation (3.71) can only be used if the damping is light and there is more than one 
overshoot. Equation (3.67) can now be employed to calculate the undamped natural 
frequency 



iO n 



2it 

~ C 2 



(3.72) 



3.7.2 Step response performance specification 

The three parameters shown in Figure 3.21 are used to specify performance in the 

time domain. 

(a) Rise time t T : The shortest time to achieve the final or steady-state value, for the 
first time. This can be 100% rise time as shown, or the time taken for example 
from 10% to 90% of the final value, thus allowing for non-overshoot response. 

(b) Overshoot : The relationship between the percentage overshoot and damping 
ratio is given in equation (3.68). For a control system an overshoot of between 
0 and 10% (1 < ( > 0.6) is generally acceptable. 

(c) Settling time t s : This is the time for the system output to settle down to within a 
tolerance band of the final value, normally between ±2 or 5%. 

Using 2% value, from Figure 3.21 

0.025 = Be-^ ts 

Invert 

50 = e Cw " t8 




Fig. 3.21 Step response performance specification. 
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Take natural logs 



giving 



In 50 = (oj n t s 




(3.73) 



The term (l/£w n ) is sometimes called the equivalent time constant T c for a second- 
order system. Note that In 50 (2% tolerance) is 3.9, and In 20 (5% tolerance) is 3.0. 
Thus the transient period for both first and second-order systems is three times the 
time constant to within a 5% tolerance band, or four times the time constant to 
within a 2% tolerance band, a useful rule-of-thumb. 



3.8 Response of higher-order systems 

Transfer function techniques can be used to calculate the time response of higher- 
order systems. 

Example 3.8 (See also Appendix 1, examp38.m) 

Figure 3.22 shows, in block diagram form, the transfer functions for a resistance 
thermometer and a valve connected together. The input x,(t) is temperature and the 
output x Q {t ) is valve position. Find an expression for the unit step response function 
when there are zero initial conditions. 

Solution 

From Figure 3.22 



Y (-) - 25 

oV> s(l + 2s)(s 2 + s + 25) 


(3.74) 


12.5 

s(s + 0.5)(.s 2 + tr + 25) 


(3.75) 


A B Cs + D 

~ S + (s + 0.5) (.s' + 0.5) 2 + (4.97) 2 


(3.76) 



Resistance Thermometer 



Valve 



=1/s 


1 

1 + 2s 




25 








s 2 + s + 25 





_ >A(s) 



Fig. 3.22 Block diagram representation of a resistance thermometer and valve. 
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Note that the second-order term in equation (3.76) has had the ‘square’ completed 
since its roots are complex ( b 2 < 4 ac). Equate equations (3.75) and (3.76) and multi- 
ply both sides by s(s + 0.5)(s 2 + s + 25). 

12.5 = (s 3 + 1.5^ 2 + 25.5s + 12.5)A + (s 3 + s 2 + 25s)B 

(3 77) 

+ (s 3 + 0.5^ 2 )C + ( s 2 + 0.5 s)D 

Equating coefficients 

(s 3 ) : 0 = A + B + C 

(s 2 ) : 0=\.5A + B + 0.5C + D 

(.s') : 0 = 25.5^ + 255 + 0.57) 

(.S’°) : 12.5 = 12.5+ 



Solving the four simultaneous equations 

A = l, B = -1.01, C = 0.01, D = —0.5 



Substituting back into equation (3.76) gives 

1 1.01 0.0b -0.5 

o(,S “ s ~ (s + 0.5) + ( S + 0.5) 2 + (4.97) 2 

Inverse transform 

x 0 (t)= 1 - 1.01e~°' 5r — 0.01e~°' 5 '(10.16sin4.97t — cos4.97t) 



(3.78) 

(3.79) 



Equation (3.79) shows that the third-order transient response contains both first- 
order and second-order elements whose time constants and equivalent time constants 
are 2 seconds, i.e. a transient period of about 8 seconds. The second-order element 
has a predominate negative sine term, and a damped natural frequency of 4.97 rad/s. 
The time response is shown in Figure 3.23. 




Fig. 3.23 Time response of third-order system. 
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3.9 Further problems 



Example 3.9 

A ship has a mass m and a resistance C times the forward velocity u(t). If the thrust 
from the propeller is K times its angular velocity ui(t), determine: 

(a) The first-order differential equation and hence the transfer function relating U ( 5 ) 
and w(j). 

When the vessel has the parameters: m = 18 000 x 10 3 kg, C = 150 000 Ns/m, and 
K = 96 000 Ns/rad, find, 

(b) the time constant. 

(c) an expression for the time response of the ship when there is a step change of u(t) 
from 0 to 12.5 rad/s. Assume that the vessel is initially at rest. 

(d) What is the forward velocity after 

(i) one minute 

(ii) ten minutes. 



Solution 



(a) m(du/dt ) + Cu = Ku(t) 



U 

U) 



(s) 



K/C 

1 + ( m/C)s 



(b) 120 seconds 

(c) u{t) = 8(1 - e- 0 00833 ') 

(d) (i) 3. 148 m/s 
(ii) 7.946 m/s 



Example 3.10 

(a) Determine the transfer function relating V 2 (. s) and V\ (s) for the passive electrical 
network shown in Figure 3.24. 

(b) When C = 2pF and R\ = R 2 = 1 Mfl, determine the steady-state gain K and 
time constant T. 

(c) Find an expression for the unit step response. 




v 2 (t) 



Fig. 3.24 Passive electrical network. 
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Solution 



V 2 Ri/ R\ + i?2 



V\ v ' l + (R i R 2 C/Ri+R 2 )s 



(b) 0.5 

1.0 seconds 



(c) v o (0 = 0.5(1 - e-0 



Example 3.11 

Determine the values of w n and £ and also expressions for the unit step response for 
the systems represented by the following second-order transfer functions 



(i) 2.0 

1 .0 (Critical damping) 
x Q (t ) = 1 - e~ 2 '(l + 2 1 ) 

(ii) 2.236 

1.342 (Overdamped) 
x 0 (r) = 2 — 2.5e~‘ + 0.5e -5 ' 

(iii) 1.0 

0.5 (Underdamped) 

x 0 (t ) = 1 - e _0 - 5 '( cos 0.866? + 0.577 sin 0.866t) 

Example 3.12 

A torsional spring of stiffness K, a mass of moment of inertia / and a fluid damper 
with damping coefficient C are connected together as shown in Figure 3.25. The 
angular displacement of the free end of the spring is 6((r) and the angular displace- 
ment of the mass and damper is 9 0 (t). 




Solution 



n 




K 




Fig. 3.25 Torsional system. 



/ 



u 

c 
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(a) Develop the transfer function relating 9i(s) and d 0 (T). 

(b) If the time relationship for dj(?) is given by dj(?) = At then find an expression for 
the time response of 9 0 {t). Assume zero initial conditions. What is the steady- 
state error between dj(?) and d 0 (?)? 

Solution 

(a) ! w= +'»*+ 1 

(b) 9 0 (t) = At — 0.2 + e~ 25 '(0.2 cos 9.682? - 0.361 sin 9.682?) 

0.2 radians. 

Example 3.13 

When a unity gain second-order system is subject to a unit step input, its transient 
response contains a first overshoot of 77%, occurring after 32.5 ms has elapsed. Find 

(a) the damped natural frequency 

(b) the damping ratio 

(c) the undamped natural frequency 

(d) the system transfer function 

(e) the time to settle down to within ±2% of the final value 

Solution 

(a) 96.66 rad/s 

(b) 0.083 

(c) 96.99 rad/s 



w w 0.106 x 10-V+ 1.712 x 10- 3 s+ 1 
(e) 0.486 seconds 



Example 3.14 

A system consists of a first-order element linked to a second-order element without 
interaction. The first-order element has a time constant of 5 seconds and a steady- 
state gain constant of 0.2. The second-order element has an undamped natural 
frequency of 4 rad/s, a damping ratio of 0.25 and a steady-state gain constant of 
unity. 

If a step input function of 10 units is applied to the system, find an expression for 
the time response. Assume zero initial conditions. 

Solution 

Xo (?) = 2.0 - 2.046e~ 0,2 ' + e~ r (o.046 cos v^? - 0.094 sin \/l5?) 




4 



Closed-loop control 
systems 



4.1 Closed-loop transfer function 

Any system in which the output quantity is monitored and compared with the input, 
any difference being used to actuate the system until the output equals the input is 
called a closed-loop or feedback control system. 

The elements of a closed-loop control system are represented in block diagram 
form using the transfer function approach. The general form of such a system is 
shown in Figure 4.1. 

The transfer function relating R(s) and C(s) is termed the closed-loop transfer 
function. 

From Figure 4.1 



C(s) = G(s)E(s) 


(4.1) 


B(s ) : //(.v)C(v) 


(4.2) 


E(s) = R(s) - B(s) 


(4.3) 



Substituting (4.2) and (4.3) into (4.1) 

C(s) = G(s){R(s) - H(s)C(s )} 
C(s) = G(s)R(s) - G(s)H(s)C(s) 
C0){1 + G(s)H(s)} = G(s)R(s) 



C (s)= 

R ; 1 + G(s)H(s) 



(4.4) 



The closed-loop transfer function is the forward-path transfer function divided by 
one plus the open-loop transfer function. 
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Forward Path 

► 



Summing 




Fig. 4.1 Block diagram of a closed-loop control system. R(s) = Laplace transform of reference input r(f); 
C(s) = Laplace transform of controlled output c(f); B(s) = Primary feedback signal, of value H(s)C(s); 
F(s) = Actuating or error signal, of value R(s) — B(s)', G(s) = Product of all transfer functions along the 
forward path; H(s) = Product of all transfer functions along the feedback path; G(s)H(s) = Open-loop 
transfer function; (g) = summing point symbol, used to denote algebraic summation; • = Signal take-off 
point; — » = Direction of information flow. 



4.2 Block diagram reduction 
4.2.1 Control systems with multiple loops 

A control system may have several feedback control loops. For example, with a ship 
autopilot, the rudder-angle control loop is termed the minor loop, whereas the 
heading control loop is referred to as the major loop. When analysing multiple loop 
systems, the minor loops are considered first, until the system is reduced to a single 
overall closed-loop transfer function. 

To reduce complexity, in the following examples the function of s notation (s) used 
for transfer functions is only included in the final solution. 

Example 4.1 

Find the closed-loop transfer function for the system shown in Figure 4.2. 

Solution 

In Figure 4.2, the first minor loop to be considered is G 1 H 3 . Using equation (4.4), this 
may be replaced by 



G 3 

1 + G3//3 



Gml = 



(4.5) 
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First Minor Loop 




Fig. 4.2 Multiple loop control system. 



Now G m i is multiplied by, or in cascade with Gi- Hence the combined transfer 
function is 



GiG m \ 



G2G3 
I +G 3//3 



(4.6) 



The reduced block diagram is shown in Figure 4.3. 

Following a similar process, the second minor loop G m 2 may be written 



Cm2 



G 2 G 3 
I + G3//3 
1 _i_ G2G3H2 
1 I+G3//3 



Multiplying numerator and denominator by 1 + G 3//3 



G m 2 



G 2 G 3 

1 + G 3//3 + g 2 g 3 h 2 



But G m 2 is in cascade with G\, hence 



GiG m2 



G1G2G3 

1 + G 3//3 + G 2 G 3//2 



(4.7) 



Transfer function (4.7) now becomes the complete forward-path transfer function as 
shown in Figure 4.4. 
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Second Minor Loop 




Fig. 4.3 First stage of block diagram reduction. 




The complete, or overall closed-loop transfer function can now be evaluated 

„ G1G2G3 _ 

Ws : +G 3 H } +G 2 G,H 2 

R 1 ' 1 1 G\G 2 G 2 H\ 

1 l+G3H 3 +G 2 G3H 2 



Multiplying numerator and denominator by 1 + G 3 H 3 + G 2 G 3 H 2 



C Gi(s)G 2 (.s)G 3 (s) 

R 1 + G 3 (s)H 3 (s) + G 2 (s)G 3 {s)H 2 (s) + G\(s)G 2 (s)G 3 {s)H ( (s) 



( 4 . 8 ) 
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4.2.2 Block diagram manipulation 

There are occasions when there is interaction between the control loops and, for the 
purpose of analysis, it becomes necessary to re-arrange the block diagram configur- 
ation. This can be undertaken using Block Diagram Transformation Theorems. 

Table 4.1 Block Diagram Transformation Theorems 
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Example 4.2 

Moving a summing point ahead of a block. 

Equation Equation 

Z = GX±Y Z={X±(l/G)Y}G (4.9) 

Z=GX±Y 

A complete set of Block Diagram Transformation Theorems is given in Table 4.1. 



Example 4.3 

Find the overall closed-loop transfer function for the system shown in Figure 4.6. 
Solution 

Moving the first summing point ahead of G\ , and the final take-off point beyond G 4 
gives the modified block diagram shown in Figure 4.7. The block diagram shown in 
Figure 4.7 is then reduced to the form given in Figure 4.8. The overall closed-loop 
transfer function is then 



C 

R 



(s) 



G\ G2G3G4 

(1+GiG 2 #i)(1+G 3 G4#2) 

1 1 01 (72 (*304 #3 

1 (0i0 4 )(l+GiG2tfi)(l+G 3 G 4 //2) 

Gi(s)G2(s)G3(s)G4(s) 

(1 + G i (5)G 2 (j)//,( 5))(1 + G 3 (s)G 4 (s)ff 2 (s)) + G 2 (s)G 3 (s)// 3 (s) 



(4.10) 




Fig. 4.5 Moving a summing point ahead of a block. 




Fig. 4.6 Block diagram with interaction. 
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Fig. 4.7 Modified block diagram with no interaction. 




Fig. 4.8 Reduced block diagram. 



4.3 Systems with multiple inputs 
4.3.1 Principle of superposition 

A dynamic system is linear if the Principle of Superposition can be applied. This 
states that ‘The response y(t ) of a linear system due to several inputs x\(t), 
x 2 (t ), . . . , x n (t), acting simultaneously is equal to the sum of the responses of each 
input acting alone’. 



Example 4.4 

Find the complete output for the system shown in Figure 4.9 when both inputs act 
simultaneously. 

Solution 

The block diagram shown in Figure 4.9 can be reduced and simplified to the form 
given in Figure 4.10. Putting AM.v) = 0 and replacing the summing point by +1 gives 
the block diagram shown in Figure 4.11. In Figure 4.1 1 note that C' 1 (s) is response to 
R\{s) acting alone. The closed-loop transfer function is therefore 



C 1 

R\ 



(s) 



G\G 2 

\+g 2 h 2 

i G\G 2 H\ 

1 ^ 1 +g 2 h 2 
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or 



CV) 



G\(s)G2(s)Ri(s) 

1 + G 2 (s)H 2 (s ) + G ] (s)G 2 (.s)H i (s) 



(4.11) 



Now if /?i (,y) = 0 and the summing point is replaced by —1, then the response C n (.s) 
to input R 2 (s) acting alone is given by Figure 4.12. The choice as to whether 
the summing point is replaced by +1 or —1 depends upon the sign at the summing 
point. 

Note that in Figure 4.12 there is a positive feedback loop. FLence the closed-loop 
transfer function relating R 2 (s) and C n (.sj is 

— (s) = [+GlH 2 

R' I f CiG 2 //i ) 

V 1 +G1H2 ) 



fll (S) + 





Fig. 4.9 System with multiple inputs. 




Fig. 4.10 Reduced and simplified block diagram. 
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Fig. 4.11 Block diagram for /?i(s) acting alone. 




Fig. 4.12 Block diagram for R 2 (s) acting alone. 



or 



r n f = 

y> 1 + G 2 (s)H 2 (s) + G x {s)G 2 {s)H [ (s) K ' 

It should be noticed that the denominators for equations (4.11) and (4.12) are 
identical. Using the Principle of Superposition, the complete response is given by 

C(,v) = C'U) + C"(j) (4.13) 



or 



C(s) 



(G l (s)G 2 (s))R l (s) - (Gi(j)Ch(j)gi(s))fl 2 (3) 
1 + G 2 (s)H 2 (s ) + Gi(s)G 2 (s)Hi(s) 



(4.14) 



4.4 Transfer functions for system elements 
4.4.1 DC servo-motors 

One of the most common devices for actuating a control system is the DC servo- 
motor shown in Figure 4.13, and can operate under either armature or field control. 

(a) Armature control : This arrangement is shown in schematic form in Figure 4. 14. 
Now air gap flux 4> is proportional to if, or 



$ = Km 



(4.15) 
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Fig. 4.13 Simple DC servo-motor, 
where Kf^ is the field coil constant. 

Also, torque developed T m is proportional to the product of the air gap flux and 
the armature current 



T m (t) = <t>/d am /a(0 



( 4 . 16 ) 
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Fig. 4.14 DC servo-motor under armature control. e a (t) = Armature excitation voltage; eb(f) = Backemf; 
/ a (f) = Armature current; R a = Armature resistance; L a = Armature inductance; e f = Constant field 
voltage; 4 = Constant field current; T m = Torque developed by motor; 9(t) = Shaft angular displacement; 
ui(t) = Shaft angular velocity = d(9/dt. 



where K am is the armature coil constant. 

Substituting (4.15) into (4.16) gives 

T m (t) = (K a K am i[)i a (t) (4.17) 

Since if is constant 

T m (t) = K.J.Jt) (4.18) 

where the overall armature constant K a is 

K a = K fd K am i { (4.19) 

When the armature rotates, it behaves like a generator, producing a back emf eb(0 
proportional to the shaft angular velocity 

e b (t) = K b ^ = K h co(t) (4.20) 

where K h is the back emf constant. 

The potential difference across the armature winding is therefore 

e a (t)-e b (t) = L a ^+R a i a (4.21) 

Taking Laplace transforms of equation (4.21) with zero initial conditions 

E a (s) - E b (s) = (L a s + U)l d (s) (4.22) 

Figure 4.15 combines equations (4.18), (4.20) and (4.22) in block diagram form. 
Under steady-state conditions, the torque developed by the DC servo-motor is 



T m (t) = {e a (t)-KMt)}§ 
Tv 
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or 



T m (t) = e,{t) - 0^) uj(t) (4.23) 

From equation (4.23), the relationship between T m (t), u>(t) and E a (t) under steady- 
state conditions is shown in Figure 4.16. 

(b) Field control'. This arrangement is shown in schematic form in Figure 4.13, 
with the exception that the armature current i a is held at a constant value. Equation 
(4. 1 7) may now be written as 



T m (t) = (K {A K am i a )i f (t) (4.24) 

and since i a is a constant, then 

T m (t) = Kfd(t) (4.25) 

where the overall field constant Kf is 

K { = K fd K am i a (4.26) 

In this instance, the back emf eb does not play a part in the torque equation, but it 
can produce difficulties in maintaining a constant armature current i a . 

The potential difference across the field coil is 

e f (t)= L f ^-+ R f k ( 4.27 ) 

at 

Taking Laplace transforms of equation (4.27) with zero initial conditions 

Eds) = (L f s + Rt)h(s) (4.28) 



Figure 4.17 combines equations (4.25) and (4.28) in block diagram form. 
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Fig. 4.16 Steady-state relationship between T m (t), ui(t) and e a (f) for an armature controlled DC servo- 
motor. 




Fig. 4.17 Block diagram representation of field controlled DC servo-motor. 



Under steady-state conditions, the torque developed by the DC servo-motor is 

Tm(t) = {jf) ef(t} (4 ' 29) 

The relationship between T m (t), Cf (?) and u{t) under steady-state conditions is shown 
in Figure 4.18. 



4.4.2 Linear hydraulic actuators 

FLydraulic actuators are employed in such areas as the aerospace industry because 
they possess a good power to weight ratio and have a fast response. 

Figure 4.19 shows a spool-valve controlled linear actuator. When the spool-valve 
is moved to the right, pressurized hydraulic oil flows into chamber (1) causing the 
piston to move to the left, and in so doing forces oil in chamber (2) to be expelled to 
the exhaust port. 

The following analysis will be linearized for small perturbations of the spool-valve 
and actuator. 



76 Advanced Control Engineering 




T m (t) 

(Nm) 



Increc 

i 


asing e,(f) 

i 














w(t) (rad/s) 
(b) 



Fig. 4.18 Steady-state relationship between T m (t), 6f(t) and w(f) for a field controlled DC servo-motor. 



Pe= 0 P s P e— 0 




Qleak 



Fig. 4.19 Spool-valve controlled linear actuator. 



It is assumed that: 

• the supply pressure P s is constant 

• the exhaust pressure P e is atmospheric 

• the actuator is in mid-position so that V\ = V 2 = V Q which is half the total volume 
of hydraulic fluid V t 

• the hydraulic oil is compressible 

• the piston faces have equal areas A 

• Q\ and Q 2 are the volumetric flow-rates into chamber (1) and out of chamber (2) 

• the average, or load flow-rate Ql has a value (Q \ + Q 2 )/ 2 

• P 1 and P 2 are the fluid pressures in chamber (1) and chamber (2) 

• the load pressure P l has a value (Pi — P 2 ) 
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(a) Actuator analysis'. The continuity equation for the chambers may be written 

Y j Qi n — ^ 2out = (rate of change of chamber volume) 

+ (rate of change of oil volume) (4.30) 



In equation (4.30), the rate of change of chamber volume is due to the piston 
movement, i.e. dV/d t. The rate of change of oil volume is due to compressibility 
effects, i.e.: 

Bulk Modulus of oil, (3 = Volumetric stress/Volumetric strain 



P 



d P 

-dV/V 



(4.31) 



Note that in equation (4.31) the denominator is negative since an increase in pressure 
causes a reduction in oil volume. 

Hence 



dV _ d P 

Giving, when differentiated with respect to time 



dV 

dr 



F\ dP 
f3J dr 



(4.32) 



For chamber ( 1 ), equation (4.30) may be expressed as 



Q\ - Sleak = — 



and for chamber (2) 



Now 



Gleak ' Ql ~ 



VA d/h 

p) dr 



dF 2 /V 2 \ dP 2 
dT + \J3 ) ~di 



(4.33) 



(4.34) 



Ql = 



Qi + Qi 



Thus, from eqations (4.33) and (4.34) 

If leakage flow-rate Q\eak is laminar, then 



- V 2 



d P : 
dr 



(4.35) 



Q\eak — C p P L 



(4.36) 
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where C p is the leakage coefficient. Also, if V\ = V 2 = V 0 , then 

dV 2 _ dVi _ dV Q 

dr dr dr 



Hence equation (4.35) can be written 



Q L 



CvPl + 



dV 0 
d t 



2(3 dr 1 



Pi) 



or 



Gl 



C p P l + A 



dAo 

dr 



V t d P L 
4(3 dr 



where 



dV ° _ A dx ° 

dr dr 



(4.37) 



(4.38) 



(4.39) 



and 



V 0 = 



Vt 

2 



(b) Linearized spool-valve analysis : Assume that the spool-valve ports are rectan- 
gular in form, and have area 

A v = WX V (4.40) 

where W is the width of the port. 

From orifice theory 



Gi 



Cd wx v J-(P S - P\ ) 



(4.41) 



and 



G2 



C<| WXy 




(Pi - 0) 



whereo Cd is a coefficient of discharge and p is the fluid density. 
Equating (4.41) and (4.42) 



(4.42) 



Ps - Pi = Pi 



(4.43) 



since 



Pl = Pi~ Pi 



Equation (4.43) may be re-arranged to give 



Pi = 



P s + Pl 



2 



(4.44) 
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From equations (4.41) and (4.42) the load flow-rate may be written as 



Ql = C d fFX v y^ ^ 2 j (4.45) 

Hence 



Q l = F(X v ,P l ) 



(4.46) 



Equation (4.45) can be linearized using the technique described in section 2.7.1. If q l, 
x v and pl are small perturbations of parameters Q l, Xy and Pl about some operat- 
ing point ‘a’, then from equation (4.46) 



</r 



<9z v 



Xyj + 
a 



gg L 

ap L 



P L 

a 



(4.47) 



or 



</L = ^qA'v - K c Pl 



(4.48) 



where 



Z q (flow gain) 



C d IF 




(P s - ^La) 



(4.49) 



and 

(flow-pressure coefficient which has a negative value) 

-c d iFx va / r 

C 2 Y p(P s - P La ) 



(4.50) 



Note that Pi_ a and Z va are the values of Pl and X v at the operating point ‘a’. 

The relationship between Q l, Pl and X v , from equation (4.45), together with the 
linearized relations q L , Pl and ,v v are shown in Figure 4.20. 

Equation (4.39) is true for both large and small perturbations, and so can be 
written 



. _ A d*o , ~ , V t d;?L 

9L - A ~dF +CpPL + 40~dF 



Equating (4.48) and (4.51) gives 



(4.51) 



KqXy 



= A d ” + (Cp + K c )p L + 



Vt dp L 
4/3 d t 



(4.52) 



Taking Laplace transforms (zero initial conditions), but retaining the lower-case 
small perturbation notation gives 



Xq A\(.V) 



F t 

Asx 0 (s) + (Cp + K c ) + —s 



Pl(s) 



(4.53) 
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Fig. 4.20 Pressure-Flow-rate characteristics fora spooi-valve. 



The force to accelerate the mass m is shown in Figure 4.21. From Figure 4.21 



^2 F x = mx 0 
AP^ = mx 0 



(4.54) 



Take Laplace transforms with zero initial conditions and using lower-case notation 



rn 2 

Pl(s) = —rx 0 (s) 



Inserting equation (4.55) into (4.53) gives 

K q x v {s) = Asx 0 (s) + |(C P + K c ) + ^.s'j (^s 2 ^jx 0 (s) 



(4.55) 



(4.56) 



Equation (4.56) may be re-arranged to give the transfer function relating x 0 (s) and 
Xy(s) 



ap l 




A 


m 



Xo (f), X 0 (f), Xo(t) 



Fig. 4.21 Free-body diagram of load on hydraulic actuator. 
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* w = a 

■ r " *{(@y+( a s s )*+ 1 } 



Equation (4.57) can be written in the standard form 

K h 



0 / \ 

— ( s ) = ~r . \ 

Xv s(-r.s' 2 + — s+ 1 ) 
V“4 “” h / 



where 



K h (hydraulic gain) 
uj nh (hydraulic natural frequency) 

Ch (hydraulic damping ratio) 




/Cp + AA / 4 p 

v 2 ) \j mViA 2 



(4.57) 



(4.58) 



Since the Bulk Modulus of hydraulic oil is in the order of 1 .4 GPa, if m and V t are 
small, a large hydraulic natural frequency is possible, resulting in a rapid response. 
Note that the hydraulic damping ratio is governed by Cp and K l: . T o control the level 
of damping, it is sometimes necessary to drill small holes through the piston. 



4.5 Controllers for closed-loop systems 
4.5.1 The generalized control problem 

A generalized closed-loop control system is shown in Figure 4.22. The control 
problem can be stated as: ‘The control action u(t) will be such that the controlled 
output c(t ) will be equal to the reference input r\(t) for all values of time, irrespective 
of the value of the disturbance input r 2 ( 0 ’- 

In practice, there will always be transient errors, but the transient period should be 
kept as small as possible. It is usually possible to design the controller so that steady- 
state errors are minimized, or ideally, eliminated. 




Fig. 4.22 Generalized closed-loop control system. 
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4.5.2 Proportional control 

In this case, the control action, or signal is proportional to the error e{t) 

u(t) = K\e{t) (4.59) 

where K\ is the proportional gain constant. 

If the plant dynamics are first-order, then Figure 4.22 can be described as shown in 
Figure 4.23. The plant transfer function is 

(U(s) - R 2 (s)) = C(s) (4.60) 

And the proportional control law, from equation (4.59) becomes 

U(s) = K l (R l (s)-C(s)) (4.61) 



Inserting equation (4.61) into equation (4.60) gives 

(1 + Ts) 

which can be written as 



C(s) = 



(4.62) 



{(1 +K t K) + Ts}C(s ) = K[KR[(s) - KR 2 (s ) 
Re-arranging equation (4.63) gives 



C(s) = 



I 1 + (i+J,/c) • s ’} 



(4.63) 



(4.64) 



When r\(t) is a unit step, and r 2 (t) is zero, the final value theorem (equation (3.10)) 
gives the steady-state response 



c(t) = 



K\K 
1 + KiK 



as t — » oo. 



R 2 (s) 



Proportional 

Controller 



«i(s) +. 




Fig. 4.23 Proportional control of a first-order plant. 
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When r 2 (?) is a unit step, and r\{t) is zero, the final value theorem (equation (3.10)) 
gives the steady-state response 

c(,)= -(rrbf) as '^°°- 

Hence, for the system to have zero steady-state error, the terms in equation (4.64) 
should be 



KiK \ 

i + k { k) 

\1 + K\KJ 



= 1 
= 0 



(4.65) 



This can only happen if the open-loop gain constant K\K is infinite. In practice this is 
not possible and therefore the proportional control system proposed in Figure 4.23 
will always produce steady-state errors. These can be minimized by keeping the open- 
loop gain constant K\K as high as possible. 

Since the closed-loop time-constant form equation (4.64) is 

<466) 

Then maintaining K\ K at a high value will reduce the closed-loop time constant and 
therefore improve the system transient response. 

This is illustrated in Figure 4.24 which shows a step change in r\ (t) followed by a 
step change in /-?(/). 



Summary 

For a first-order plant, proportional control will always produce steady-state errors. 
This is discussed in more detail in Chapter 6 under ‘system type classification’ where 
equations (6.63)— (6.65) define a set of error coefficients. Increasing the open-loop 




Fig. 4.24 Step response of a first-order plant using proportional control. 
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gain constant (which is usually achieved by increasing the controller gain K\) will 
reduce, but not eliminate them. A high controller gain will also reduce the transient 
period. However, as will be shown in Chapters 5 and 6, high open-loop gain 
constants can result in the instability of higher-order plant transfer functions. 



4.5.3 Proportional plus Integral (PI) control 



Including a term that is a function of the integral of the error can, with the type of 
plant shown in Figure 4.23, eliminate steady-state errors. 

Consider a control law of the form 



u(t) = K\e{t) + K 2 



edt 



Taking Laplace transforms 



U(s) 




(4.67) 



= A '( l+ i|) £(s) 

= * 1 ( I + E(s) (4.68) 

In equation (4.68), 7) is called the integral action time, and is formally defined as: 
‘The time interval in which the part of the control signal due to integral action 
increases by an amount equal to the part of the control signal due to proportional 
action when the error is unchanging’. (BS 1523). 

Inserting the PI control law given in equation (4.68) into the first-order plant 
transfer function shown in equation (4.60) gives 



C(s) 



(Ki( 1 + l/Tis)(Ri(s) - C(s)) - R 2 (s))K 

(1 + 7, s) 



(4.69) 



which can be written as 



{7,7s 2 + 7,(1 + K X K) S + K\K}C{s) = K X K( 1 + 7„s)/?i(,s) - K x KT iS R 2 (s) (4.70) 
Re-arranging gives 



C{s) 



(1 + Tjs)R\(s) — T,sR 2 (s) 

(Hi) 5,2 + T > ( ! + ih) s + 1 



(4.71) 



The denominator is now in the standard second-order system form of equation 
(3.42). The steady-state response may be obtained using the final value theorem 
given in equation (3.10). 



c(t ) = (1 + 0)ri(t) — ( 0 )r 2 (t) as t — > 00 



(4.72) 
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Fig. 4.25 Step response of a first-order plant using PI control. 

When there are step changes in r\ (t) and rjit): 

C(,) = ( ' +0W|,| - ( ||*M 
£ £ 

= Ri(s) 

c{i) = n{f) (4.73) 

Thus, when r\{t) and ri(l) are unchanging, or have step changes, there are no steady- 
state errors as can be seen in Figure 4.25. The second-order dynamics of the closed- 
loop system depend upon the values of 7), T , K\ and K. Again, a high value of K\ will 
provide a fast transient response since it increases the undamped natural frequency, 
but with higher order plant transfer functions can give rise to instability. 

Summary 

For a first-order plant, PI control will produce a second-order response. There will be 
zero steady-state errors if the reference and disturbance inputs r\(t) and r 2 {t) are 
either unchanging or have step changes. The process of including an integrator within 
the control loop to reduce or eliminate steady-state errors is discussed in more detail 
in Chapter 6 under ‘system type classification’. 

Example 4.5 (See also Appendix 1, examp45.m) 

A liquid-level process control system is shown in Figure 4.26. The system parameters 
are 

A = 2 m 2 Rf = 1 5 s/m 2 

H\ = 1 V/m K v = 0.1 ni 3 /sV K\ = 1 (controller again) 

(a) What are the values of T\ and ( when the undamped natural frequency w n is 
0.1 rad/s? 

(b) Find an expression for the time response of the system when there is a step 
change of h&(t) from 0 to 4 m. Assume zero initial conditions. 
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Fig. 4.26 Liquid-level process control system. 



The controller is given in equation (4.68). The inflow to the tank is 

Vi (t) = K v u(t) (4.74) 

The tank dynamics are expressed, using equation (2.63) as 

vi(t) - v 2 (t) = (4.75) 

and the linearized outflow is 



V2 (0 = 



Kit) 

Rt 



(4.76) 



The measured head h m (t) is obtained from the pressure transducer 

h m (t) = Hih a (t) (4.77) 

From equations (4.75) and (4.76), the tank and outflow valve transfer function is 

R[ 



Vi W 1 + AR f s 



(4.78) 



The block diagram for the control system is shown in Figure 4.27. From the block 
diagram, the forward-path transfer function G(s ) is 



G(s) = 



KiK v 



*(l+5b) 

(1 + AR[S) 
K { K,Rf(\ + T t s) 



(4.79) 



T,y(l + AR r ) 
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PI Controller Tank and valve 




Fig. 4.27 Block diagram for liquid-level process control system. 

Using equation (4.4), the closed-loop transfer function becomes 



KiK v Rf(l+TiS) 

(AR[TiS 2 +Tis) 



h * } i+ g wsssr 



(4.80) 



which simplifies to 



KiK.Rdl + T-,s) 



^ (s) = 

H d v (AR f T\)s 2 + 7,(1 + K\ K v R ( rh > + K { K,R<II\ 



(4.81) 



Equation (4.81) can be expressed in the standard form of equation (3.42) for a 
second-order system. 

Putting H\ = 1, then 



H * f , 

#d * (AIL 



(1 + T iS ) 



(a, aQ ^ (a, A' v Af + l) S + ^ 



(4.82) 



(a) Comparing the denominator terms with the standard form given in equation (3.42) 



AT i 
KiK v 

1 



1 






From equation (4.83) 



T\ = 



K t K v R f 
KiK v 1x0.1 



+ 1 = 



2C 



LulA 0.1 2 x 2 



UJ n 



= 5 seconds 



(4.83) 

(4.84) 



From equation (4.84) 



c = 



uj n T { / 1 

2 \K\K W R[ 



+ 1 



0.1 x 5 



1 



1 x 0.1 x 15 



+ 1 = 0.417 
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(b) Inserting values into equation (4.82) 

H U,\- (1 + 5 *) 

H d ’ lOO.v 2 + 8.34.9 + 1 



(4.85) 



For a step input of height 4 m 



4 

5 

Expanding by partial fractions using 3.2.4 (iv) 



H a (s) = 



0.01(1 + 55) 



0.0834.V + 0.01 



HJs) 



0.04 + 0.25 
5(5 2 + 0.0834.9 + 0.01) 



A Bs+C 

7 + 5 2 + 0.08345 + 0.01 



Multiplying through by s(s 2 + 0.08345 + 0.01) 



0.04 + 0.25 = A(s 2 + 0.08345 + 0.01) + Bs 2 + Cs 



(4.86) 



Equating coefficients 



giving 



(5 2 ) : 0 = A + B 

(5 1 ) : 0.2 = 0.0834+ + C 

(5°) : 0.04 = 0.01+ 

+ = 4 B = - 4 C = -0.1336 



Substituting values back into (4.86) and complete the square to give 



rr/x 4 -45- 0.1336 

/ f+v ) — — I d 

5 (5 + 0.041 7) 2 



0.0909 2 

Inverse transform using Laplace transform pairs (9) and (10) in Table 3.1. 



(4.87) 



H a (s) = — 



45 



5 ((5 + 0.0417) 2 + 0.0909 2 



- 1.4697 



0.0909 



1(5 + 0.041 7) 2 + 0.0909 2 



h a (t) = 4- 4e-° 0417 ' ^cos0.0909t - sin 0.0909^ 



- 1.4697e-° 0417r sin 0.0909? 



which simplifies to give 

hjt) = 4[1 - e-°' 0417i (cos0.0909t- 0.0913 sin 0.09090] (4.88) 
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Fig. 4.28 Response of the PI controlled liquid-level system shown in Figure 4.26 to a step change in /ld(f) 
from 0 to 4 m. 



In equation (4.88) the amplitude of the sine term is small, compared with the cosine 
term, and can be ignored. Hence 

/z a (0 = 4(1 - e’ 00417 ' cos 0.0909?) (4.89) 

The time response depicted by equation (4.89) is shown in Figure 4.28. 



4.5.4 Proportional plus Integral plus Derivative (PID) control 

Most commercial controllers provide full PID (also called three-term) control action. 
Including a term that is a function of the derivative of the error can, with high-order 
plants, provide a stable control solution. 

Proportional plus Integral plus Derivative control action is expressed as 



f de 

u(t) = K\e(t) + K 2 edt + K 3 — 



(4.90) 



Taking Laplace transforms 



U(s)= ( K\ + — + K^s ) E(s) 

= Kl[l+ ^s + f l s]m 

= K l ( 1 + ^ + T a s ) E(s) 



(4.91) 
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In equation (4.91), T d is called the derivative action time, and is formally defined as: 
‘The time interval in which the part of the control signal due to proportional action 
increases by an amount equal to the part of the control signal due to derivative action 
when the error is changing at a constant rate’ (BS 1523). 

Equation (4.91) can also be expressed as 



U(s) 



K\{T t T d s 2 + TjS+l) 
T,s 



E(s) 



(4.92) 



4.5.5 The Ziegler-Nichols methods for tuning PID controllers 

The selection of the PID controller parameters K\, 7) and T d can be obtained using 
the classical control system design techniques described in Chapters 5 and 6. In the 
1940s, when such tools were just being developed, Ziegler and Nichols (1942) devised 
two empirical methods for obtaining the controller parameters. These methods are 
still in use. 

(a) The Process Reaction Method: This is based on the assumption that the open- 
loop step response of most process control systems has an S-shape, called the process 
reaction curve, as shown in Figure 4.29. The process reaction curve may be approxi- 
mated to a time delay D (also called a transportation lag) and a first-order system of 
maximum tangential slope R as shown in Figure 4.29 (see also Figure 3.13). 

The Process Reaction Method assumes that the optimum response for the closed- 
loop system occurs when the ratio of successive peaks, as defined by equation (3.71), 
is 4:1. From equation (3.71) it can be seen that this occurs when the closed-loop 
damping ratio has a value of 0.21. The controller parameters, as a function of R and 
D, to produce this response, are given in Table 4.2. 




Fig. 4.29 Process reaction curve. 
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Table 4.2 Ziegler-Nichols PID parameters using the 
Process Reaction Method 



Controller type 


Ki 


T; 


T A 


p 


l/RD 


- 


— 


pi 


0.9/RD 


D/ 0.3 


- 


PID 


1.2/RD 


2D 


0.5 D 


Table 4.3 Ziegler-Nichols PID parameters using the 


Continuous Cycling Method 






Controller type 


Kx 


T, 


T A 


p 


KJ2 


— 


— 


pi 


KJ2.2 


TJX. 2 


- 


PID 


KJX.l 


TJ2 


TJ8 



Note that the Process Reaction Method cannot be used if the open-loop step 
response has an overshoot, or contains a pure integrator(s). 

(b) The Continuous Cycling Method : This is a closed-loop technique whereby, 
using proportional control only, the controller gain K\ is increased until the system 
controlled output c(t) oscillates continually at constant amplitude, like a second- 
order system with no damping. This condition is referred to as marginal stability 
and is discussed further in Chapters 5 and 6. This value of controller gain is called 
the ultimate gain K u , and the time period for one oscillation of c(t) is called the 
ultimate period T u . The controller parameters, as a function of K u and T u , to provide 
a similar closed-loop response to the Process Reaction Method, are given in 
Table 4.3. 

The two Ziegler-Nichols PID tuning methods provide a useful ‘rule of thumb’ 
empirical approach. The control system design techniques discussed in Chapters 5 
and 6 however will generally yield better design solutions. 

Of the two techniques, the Process Reaction Method is the easiest and least 
disruptive to implement. In practice, the measurement of R and D is very subjective, 
and can lead to errors. 

The Continuous Cycling Method, although more disruptive, has the potential to 
give better results. There is the risk however, particularly with high performance 
servo-mechanisms, that if K n is increased by accident to slightly above the marginal 
stability value, then full instability can occur, resulting in damage to the system. 

Actuator saturation and integral wind-up 

One of the practical problems of implementing PID control is that of actuator 
saturation and integral wind-up. Since the range of movement in say, a control valve, 
has physical limits, once it has saturated, increasing the magnitude of the control 
signal further has no effect. However, if there is a difference between desired and 
measured values, the resulting error will cause a continuing increase in the integral 
term, referred to as integral wind-up. When the error term changes its sign, the 
integral term starts to ‘unwind,’ and this can cause long time delays and possible 
instability. The solution is to limit the maximum value that the integral term 
can have. 
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4.5.6 Proportional plus Derivative (PD) control 

Proportional plus Derivative control action is expressed as 

de 

u(t) = K\ e( t ) + K 3 , (4.93) 

at 

Taking Laplace transforms 

U(s) = K^\+^j E(s) 

= K x {\ + T d s)E(s) (4.94) 

The inclusion of a derivative term in the controller generally gives improved damping 
and stability. This is discussed in more detail in Chapters 5 and 6. 



4.6 Case study examples 



Example 4.6.1 CNC Machine-Tool Positional Control (See also Appendix 1, 
examp461.m) 

The physical configuration and block diagram representation of a CNC machine- 
tool is shown in Figures 1.10 and 1.11. The fundamental control problem here is that, 
by design, the lead-screw (by the use of re-circulating ball-bearings) is friction-free. 
This means that the positional control system will have no damping, and will oscillate 
continuously at the undamped natural frequency of the closed-loop system. 

Damping can be introduced in a number of ways: 

(a) A dashpot attached to the lead-screw. This is wasteful on energy and defeats the 
objective of a friction-free system. 

(b) Velocity feedback : A signal from a sensor that is the first derivative of the 
output (i.e. velocity) will produce a damping term in the closed-loop transfer 
function. 

(c) PD control : A PD controller will also provide a damping term. Flowever, the 
practical realization will require an additional filter to remove unwanted high 
frequency noise (see Chapter 6 for further details on lead-lag compensation). 

Most machine-tool manufacturers employ velocity feedback to obtain the necessary 
damping. Since overshoot in a cutting operation usually cannot be tolerated, the 
damping coefficient for the system must be unity, or greater. 

For this study, the machine-tool configuration will be essentially the same as 
shown in Figure 1.10, with the exception that: 

(i) A gearbox will be placed between the servo-motor and the lead-screw to provide 
additional torque. 

(ii) The machine table movement will be measured by a linear displacement trans- 
ducer attached to the table. This has the advantage of bringing the table ‘within 
the control-loop’ and hence providing more accurate control. 
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System element dynamic equations : With reference to Figures 1.11 and 4.31 

1. Controller 



Proportional control, gain A)(V/m) 
Control signal LhCs) = K\{X<±{s) — W m (Y)) 



2. Power amplifier 

Gain AS(V/V) 

Control signal U 2 (s) = K 2 (Ui(s) — B 2 (s)) 



(4.96) 



3. DC servo-motor-. Field controlled, with transfer function as shown in Figure 4.17. 
It will be assumed that the field time constant Lf/Rf is small compared with the 
dynamics of the machine table, and therefore can be ignored. Hence, DC servo- 
motor gain Kt, (Nm/V). 

Motor Torque T m (s) = Ah 63 (.v) (4.97) 

4. Gearbox, lead-screw and machine-table : With reference to Figure 2.9 (free-body 
diagram of a gearbox), the motor-shaft will have zero viscous friction C m , hence 
equation (2.22), using Laplace notation, becomes 

X(s) = - a (T m (s) - I m s 2 d m (s)) (4.98) 

The output shaft in this case is the lead screw, which is assumed to have zero moment 
of inertia I Q and viscous friction C 0 . The free-body diagrams of the machine-table 
and lead-screw are shown in Figure 4.30. 

For lead-screw 

Work in = Work out 
bX(t)0 o {t) = F(t)x 0 (t) 
or 

F(t) = bX(t)^- (4.99) 

Xo(t ) 



X a (t), X a (t) 

m 







Fig. 4.30 Free-body diagrams of lead-screw and machine-table. 
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Now the pitch p of the lead-screw is 



o o 
K ^ 

II 


(4.100) 


Substituting (4.100) into (4.99) 




II 


(4.101) 


p 




The equation of motion for the machine-table is 




F(t) = mx o 


(4.102) 


Equating (4.101) and (4.102) gives 




X(t) = ^ (pmx 0 ) 




Taking Laplace transforms 




X(s) = ^( pms 2 X 0 (s )) 


(4.103) 


Equating (4.98) and (4.103) gives 




pms 2 X 0 (s) = ^ ( T m (s ) - I m s 2 6 m (s)) 


(4.104) 


Now 




b/a = gear ratio n 
0m(s) = n0 o (s) 




Hence 




s 2 0 m (s) = ns 2 0 o (s) 




and 




, X 0 (s) 

0 o (s) = 


(4.105) 



Equation (4.105) can be substituted into (4.104) 

7 n 7 

pms X 0 (s) = nT m (s ) - nl m -s X Q (s) 

P 



or 



nT m (s) 




(4.106) 



giving the transfer function for the gearbox, lead-screw and machine-table as 



X 0 

T m 



(s) 



n 

(pm + n 2 I m /p)s 2 



(4.107) 



where the term n 2 I m / P may be considered to be equivalent mass of 7 m referred to the 
machine-table. 
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5. Tachogenerator 



Gain H 2 (Vs/rad) 
Feedback signal B 2 (s) = H 2 s0 o (s ) 



or, from equation (4.105) 



Ih(s) 



— sX 0 {s) 
P 



6. Position transducer 

Gain H\ (V/m) 

Feedback signal X m (s) = H\X 0 (s) 



(4.108) 



(4.109) 



(4.110) 



The system element dynamic equations can now be combined in the block diagram 
shown in Figure 4.31. Using equation (4.4), the inner-loop transfer function is 



G(s) 



K 2 K 2 np 

(p 2 m + n 2 I m )s + K 2 K 2 nH 2 



(4.111) 



Again, using equation (4.4), the overall closed-loop transfer function becomes 



Ao K\K 2 K 2 np 

Xd S ( p 2 m + n 2 I m )s 2 + K 2 K 2 nH 2 s + K\K 2 K 2 npH\ 



(4.112) 



which can be written in standard form 



To 

Xd 



(s) 






p 2 m+n 2 I m 

KiK 2 K 3 npH, 




V+ 1 



(4.113) 




Position Transducer 



Fig. 4.31 Block diagram of CNC machine-tool control system. 
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Specification : The CNC machine-table control system is to be critically damped with 
a settling time of 0.1 seconds. 

Control problem : To select the controller gain K\ to achieve the settling time and 
tachogenerator constant to provide critical damping. 

System parameters 

K 2 = 2 V/V Kj = 4Nm/V 
n = 10 : 1 p = 5 x 10~ 3 m 

m = 50 kg I m = 10 x 10~ 6 kgm 2 

Hi = 60 V/m 

Calculation of K\ : In general, the settling time of a system with critical damping is 
equal to the periodic time of the undamped system, as can be seen in Figure 3.19. 
This can be demonstrated using equation (3.62) for critical damping 

x 0 (t) = [1 - e~“"'(l + uV)] 

when 



t = 27r/w n 

x 0 {t) = [1 - e- 2 "(l + 2 tt)] 

= 0.986 (4.114) 



Thus, for a settling time of 0.1 seconds for a system that is critically damped, the 
undamped natural frequency is 



LO n 



— = 62.84 rad/s 
0.1 



(4.115) 



Comparing the closed-loop transfer function given in equation (4.113) with the 
standard form given in (3.42) 



f KffK^npHC 
V P 2,n + ” 2/ m / 



(4.116) 



Hence 



_ (jfm + n 2 I m )u>l 
1 “ K 2 K 2 npH\ 

= {(5 x 10~ 3 ) 2 x 50} + (10 2 x 10 x 1 0 6 ) 1 2 

(2 x 4 x 10 x 5 x 10 -3 x 60) 

= 0.365 V/V (4.117) 

Again, comparing equation (4. 113) with the standard form (3.42) 

2C H 2 



Wn KipHi 



(4.118) 
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Hence 

H 2(K { pH | 

w n 

2 x 1 x 0.365 x 5 x 10~ 3 x 60 
= 62.84 

= 3.485 x 10~ 3 Vs/rad (4.119) 



Example 4.6.2 Temperature control system (See also Appendix 1, examp462.m) 
The general form of a temperature control system is shown in Figure 1.6 with the 
corresponding block diagram given in Figure 1.7. 

The system variables are: 

8 d (t) = Desired temperature (°C) 

6 m ( t) = Measured temperature (V) 

9 0 (t) = Actual temperature (°C) 

9 s (t) = Temperature of surroundings (°C) 
u( t) = Control signal (v) 
i>( t) = Gas flow-rate (m 3 /s) 

Qi(t) = Heat flow into room (J/s = W) 

Q 0 (t ) = Heat flow though walls (W) 

System equations 

1. Controller. The control action is PID of the form given in equation (4.91) 

U(s) = K^ 1 + ^ + T d .?j (6 d (s) - 9 m {s)) (4. 120) 

2. Gas solenoid valve : This is assumed to have first-order dynamics of the form 



V 

~U 



(s) 



K 2 

1 + T\s 



where K 2 is the valve constant (m 3 /s V). 



3. Gas burner. This converts gas flow-rate v(t) into heat flow Qi(t) i.e.: 



(4.121) 



Q,(s) = K 3 V (s) (4.122) 

where K 2 is the burner constant (Ws/m 3 ). 

4. Room dynamics : The thermal dynamics of the room are 

c\0 

Q,U) - Qo(t) = Ct^ (4.123) 

Equation (4.123) is similar to equation (2.54), where Cj is the thermal capacitance 
of the air in the room. 
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The heat flow through the walls of the building is as given in equation (2.53), i.e. 

(0o (0 - 0.(0) 



Go(0 = 



Ri 



where Rj is the thermal resistance of the walls, see equation (2.47). 
Substituting equation (4.124) into (4.123) gives 

'0o(O - 0s(O N 



Gt(0 - 

Multiplying through by Rj 



Ri 



r Ad 0 
= Cj ~di 



RtQM + 0.(0 = 00 (0 + RjCj 



ddo 

At 



(4.124) 



(4.125) 



(4.126) 



Taking Laplace transforms 

RjQAs) + 0,(0 = (1 + R t C t s)0 o (s) (4.127) 

Equation (4. 127) can be represented in block diagram form as shown in Figure 4.32. 
5. Thermometer. The thermometer equation is 

9 m {s) = H x e o { S ) (4.128) 

The complete block diagram of the temperature control system is shown in Figure 
4.33. 

From Figure 4.33 



K^KiRjjT^s 2 + TjS + l)(0 d (Q - H x 0 O (Q) 
Ts(\ + T lS ) 



+ 0,(0 = (1 + R t C t s)0 o (s) (4. 129) 



e s (s) 



Oi(s) 






L+ 


1 


0o (s) 








1 + RjCjS 





Fig. 4.32 Block diagram representation of the thermal dynamics of the room. 




Fig. 4.33 Block diagram of temperature control system. 
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Equation (4.129) can be re-arranged to give 

ir, (7) T d s 2 + T lS + l)0 d (s) + m 



0 o (s) = 






wm + i)^+i 



where the forward-path gain Kp is 

K f = K t K 2 K 3 R T 



(4.130) 



(4.131) 



Control problem-. Given the system parameters, the control problem is to determine 
the controller settings for K \ , T\ and T d . This will be undertaken using the Zeigler- 
Nichols process reaction method described in Section 4.5.5(a). 

System parameters: 

K 2 Kt, = 5 W/V Rj = 0.1 Ks/J 
Cj = 80 J/K Hi = 1.0 V/K 
T\ = 4 seconds 



Process reaction curve : This can be obtained from the forward-path transfer function 



0 o ( _ K 2 K 3 R t 
U {S) (1 + 7,^(1 + RjCjs) 



(4.132) 



Inserting values into equation (4.132) gives 



0 o ( , 

c7 W = (T 



0.5 



4.0(1 + 8 s) 



(4.133) 



Figure 4.34 shows the response to a unit step, or the process reaction curve. 

From the R and D values obtained from the process reaction curve, using the 
Zeigler-Nichols PID controller settings given in Table 4.2 



K\ = 1.2 /RD = 26.144 
7) = 2D = 3.0 seconds 
7~d = 0.5 D = 0.75 seconds 



Assuming that the temperature of the surroundings 9 s (t) remains constant, the 
closed-loop transfer function (using equation (4.130)) for the temperature control 
system, is 

— (s) = (2 - 25 * 2 + 1 * +1) (4 134) 

6 d K ’ 7.344.? 3 + 5.004s 2 + 3.229s +1 ' ’ 

The response to a step change in the desired temperature of 0-20 °C for the closed- 
loop transfer function given by equation (4.134) is shown in Figure 4.35. 

From Figure 4.35, the ratio of successive peaks is 



Cl 2 



8.92 

F98 



= 4.5 



(4.135) 
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Fig. 4.34 Process reaction curve for the temperature control system shown in Figure 4.33. 




Fig. 4.35 Closed-loop step response of temperature control system using PID controller tuned using 
Zeigler- Nichols process reaction method. 

This corresponds to a damping ratio of 0.23. These values are very close to the 
Zeigler-Nichols optimum values of 4.0 and 0.21 respectively. 

Example 4.6.3 Ship Autopilot (See also Appendix 1, examp463.m ) 

A ship has six degrees-of-freedom, i.e. it is free to move in six directions simultan- 
eously, namely three linear motions - surge (forward), sway (lateral) and heave 
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(vertical) together with three rotational motions - roll, pitch and yaw. This analysis 
considers rotation about the yaw axis (i.e. heading control) only. 

Figure 1.12 shows a typical ship autopilot system and Figure 1.13 shows the 
corresponding block diagram. Rotation about the yaw axis is in effect rotation about 
the z, or vertical axis of the vessel, called the V’ direction since r is the symbol for 
yaw-rate. Hence hydrodynamic coefficients for the yaw axis are therefore given the 
subscript V. Yaw hydrodynamic coefficients are given the symbol ' N\ In this 
analysis the dynamics of the steering gear are neglected. 

The system variables are 

\E' d (?) = Desired heading (radians) 
v]> a (?) = Actual heading (radians) 

S( t) = Rudder angle (radians) 

Figure 4.36 shows the hull free-body diagram. Disturbance effects (wind, waves and 
current) are not included. 

System equations 

1. Hull dynamics’. In Figure 4.36, X 0 Y 0 is the earth co-ordinate system where X Q is 
aligned to north. All angles are measured with respect to X Q . A consistent right- 
hand system of co-ordinates is used, with the exception of the rudder-angle, which 
has been selected to be left-hand to avoid negative coefficients in the hull transfer 
function. 




Fig. 4.36 Free-body diagram of ship hull dynamics. 
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From Figure 4.36 the equation of motion for the yaw axis is 

X> G = / Z &(0 

(«)) ~ (N r ip a (t)) = / z t/> a (0 



Taking Laplace transforms 

Ng6(s) = (I z s 2 + N r s)ip a (s) 
Flence the hull transfer function becomes 



Ip a 
5 



co 



n 6 

s(I z s + N r ) 



(4.136) 



(4.137) 



(4.138) 



2. Control action-. In this case, the autopilot (controller) is considered to provide 
proportional control only. 



8(s) = Ki(tpi(s) - t/>aC0) 



(4.139) 



3. Gyro-compass'. This provides a measured heading proportional to the actual 
heading 



ipm(s) = HvCAs) 



(4.140) 



Combining equations (4.138), (4.139) and (4.140) produces the block diagram shown 
in Figure 4.37. 

Using equation (4.4), the closed-loop transfer function is 



K,K 2 N s 

s(I z s+N t ) 

1 , K 1 K 1 N s H i 
1 + s(I z s+N r ) 



Equation (4.141) simplifies to give 



Hi 



V’d ( h 'l s 2 i ( N± A , i 

\KiK 1 N l ,Hi) !> ^ \KiK 2 N s Hi)^ ^ 1 



(4.141) 



(4.142) 




Fig. 4.37 Block diagram of ship autopilot control system. 
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Equation (4.142) is in the standard form given in equation (3.42). 

Control problem : For a specific hull, the control problem is to determine the autopilot 
setting ( K \ ) to provide a satisfactory transient response. In this case, this will be when 
the damping ratio has a value of 0.5. Also to be determined are the rise time, settling 
time and percentage overshoot. 

System parameters'. The ship to be controlled is a cargo vessel of length 161 m with a 
MARINER type hull of total displacement 17 000 tonnes. 

K 2 = 1.0 rad/V N s = 80 x 10 6 Nm/rad 

N r = 2 x 10 9 Nms/rad I z = 20 x 10 9 kgm 2 
Hi = l.OV/rad 



Inserting values into equation (4.142) gives 






(s) 



which simplifies to 




(4.143) 



(4.144) 



Comparing equation (4.144) with the standard form given in equation (3.42) 



1 250 

u>l Ki 



(4.145) 



2C_ 25 
W n Ki 



(4.146) 



Given that ( = 0.5, then from equation (4.146) 

Kx 
A 

Substituting (4.147) into (4.145) gives 

25 2 

Kf 



25 

250 



(4.147) 



Hence 



K\ = 



625 

250 



= 2.5 



(4.148) 



and from equation (4.147), 
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Fig. 4.38 Unit step response of ship autopilot control system. RiseTime (to 95%) = 23 seconds; 
Percentage Overshoot = 16.3%; 

Settling time (to ±2%) = 81 seconds. 



From equations (3.58) and (3.59) the unit step response for the ship autopilot control 
system is given by the expression 

V> a (0= 1 -e~ 0 05r ( cos 0.0866f + 0.577 sin 0.0866?) (4.150) 

Figure 4.38 shows the system unit step response. From Figure 4.38 



4.7 Further problems 



Example 4.7 

For the block diagrams shown in Figure 4.39 find an expression for the complete 
output when all inputs act simultaneously. 

Solutions 

, \ rl n _ (gj (s)G 2 (s)Gj(s))Ri (s) + (? 3 CgK 1 + G 2 (s)H 3 (s))R 2 (s) 

(a; + G 3 (s)H 2 (s) + G 2 (s)H 3 (s) + G l (s)G 2 (s)G 3 (s)H l (s) 

,, v r( . _ (gj (tj)G 2 U)G3(s)G*(s))Ri (s) - (Cl MM7 - (G 3 (s)G 4 {s))Ri(s) 

l + G 3 (s)H 2 (s) + Gi(s)G2{s)(hQ!)G4(s)Hi(,s) 

Example 4.8 

The speed control system shown in Figure 4.40 consists of an amplifier, a field- 
controlled DC servomotor and a tachogenerator. The load consists of a mass of 
moment of inertia I and a fluid damper C. The system parameters are: 

/ = 0.75 kg m 2 C = 0.5Nms/rad 
K 2 = 5 Nm/A Hi = 0.1 Vs/rad 
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(a) 




Fig. 4.39 Block diagrams for closed-loop systems. 

Amplifier l a = constant 




Find the value of K\ to give the system a closed-loop time constant of one second. 
What is the steady-state value of uj Q (t) when v;(f) has a value of 10 V. 



Solution 
0.5 A/V 
33.3 rad/s 



to 
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Fig. 4.41 Angular positional control system. ki = Error detector gain (V/rad); K 2 = Amplifier gain (A/V); 
K3 = Motor constant (Nm/A); n = Gear ratio; Hi = Tachogenerator constant (V s/rad); / G = Load moment 
of inertia (kg m 2 ); C 0 = Load damping coefficient (Nms/rad). 



Example 4.9 

Find an expression for the closed-loop transfer function relating 6\{s ) and 0 o (s ) for 
the angular positional control system shown in Figure 4.41. 



Solution 



7 rO) = 






( C a +K 2 K^H t \ 
\ k k,K.n I 



s+ 1 



Example 4.10 

A hydraulic servomechanism consists of the spool-valve/actuator arrangement 
shown in Figure 4.19 together with a ‘walking beam’ feedback linkage as shown in 
Figure 4.42. The spool-valve displacement x v (t) is given by the relationship 




Fig. 4.42 Hydraulic servomechanism with 'walking beam' feedback linkage. 
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Fig. 4.43 Block diagram of a servomechanism. 



- (fTT t) xM 

If the forward-path transfer function is given by equation (4.57), find an expression for 
the closed-loop transfer function relating X { (s) and A 0 (s). The system parameters are 

/n = 50 kg F t = 4xl0~ 3 m 3 /3=1.4GPa 

A = 0.01 m 2 K q = 10.0 m 2 /s K c = 6 x 10~ 9 m 5 /Ns 

C p = 6 x 10~ 9 m 5 /Ns a = b = 0.15m 



Solution 



i'o,, 

* (i) 



500 

0.357 x 10- 6 s 3 + 0. 12 x lO- 3 * 2 + s + 500 



Example 4.11 

A servomechanism consists of an amplifier, actuator and sensor as shown in block 
diagram form in Figure 4.43. If the input to the system is a constant velocity of the 
form 



0i(O = 2 1 

find an expression for the time response of the system. 
Solution 



f) 0 (t) = 0.667 1 - 0.0222 + e"'(0.0222 cos 7.68/ - 0.083 sin 7.68/) 



Example 4.12 

Figure 4.44 shows the elements of a closed-loop temperature control system. A 
proportional controller compares the desired value 9[(t) with the measured value 
v 0 (?) and provides a control signal u(t ) of K\ times their difference to actuate the valve 
and burner unit. The heat input to the oven Qi{t) is A3 times the control signal. 
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Proportional 




Fig. 4.44 Closed-loop temperature control system. 



The walls of the oven have a thermal resistance Rj and the oven has a thermal 
capacitance Cj and operating temperature 9 0 (t). The heat transfer equation for the 
oven may be written 



QM 



0 o (t) _ dd a 
R t T dr 



The thermometer and measurement system feed a measured value of H \ times 9 0 (t) to 
the controller. The system parameters are 



Ki = 5 K 2 = 1.5 J/V Hi = 1 V/K 
Rj = 2 K/J C T = 25 Js/W 



Find 

(a) The open-loop time constant 

(b) The closed-loop time constant 

(c) The percentage steady-state error in the output when the desired value is con- 
stant. 

Solution 

(a) 50 seconds 

(b) 3.125 seconds 

(c) 6.25% 

Example 4.13 

Figure 4.45 shows the block diagram representation of a process plant being con- 
trolled by a P1D controller. 

(a) Find an expression for the complete response C(.v) when R \ (.v) and R 2 (s) act 
simultaneously. 

(b) Using the Ziegler-Nichols Process Reaction Method, determine values for K\, T, 
and 7d when T\ = 10 seconds and T 2 = 20 seconds. 
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Fig. 4.45 Process plant under PID control. 



(c) Insert the values into the expression found in (a). Using MATLAB, or otherwise, 
determine the response c(t) when r\{t) is a unit step and r 2 (t) is zero. What is the 
ratio of successive peaks? 

Solution 



(a) C(s) 



KiiT.TiS 2 + T,s + l)tf,(s) + TisR 2 (s) 

Tih T 2 s 3 + T,(T\ + T 2 + K ] T d ) s 2 + T,(l + Ki)s + K\ 



(b) K\ = 17.2 T\ = 6 seconds T d = 1.5 seconds 



(c) C(s) 



ll.2(9s 2 + 6s+l) 

^(1200U + 334.8s 2 + 109.2^+ 17.2) 



2.75 



5 



Classical design in the 
s-plane 



5.1 Stability of dynamic systems 

The response of a linear system to a stimulus has two components: 

(a) steady-state terms which are directly related to the input 

(b) transient terms which are either exponential, or oscillatory with an envelope of 
exponential form. 

If the exponential terms decay as time increases, then the system is said to be stable. If 
the exponential terms increase with increasing time, the system is considered unstable. 
Examples of stable and unstable systems are shown in Figure 5.1. The motions 
shown in Figure 5.1 are given graphically in Figure 5.2. (Note that (b) in Figure 
5.2 does not represent (b) in Figure 5.1.) The time responses shown in Figure 5.2 can 
be expressed mathematically as: 

For (a) (Stable) 





x a (t) = Ae al sin(u;? + </>) 


(5.1) 


For (b) (Unstable) 


x Q (t) = Ae at sin (cet + </>) 


(5.2) 


For (c) (Stable) 


b 

1 

<U 

II 

O 

K 


(5.3) 


For (d) (Unstable) 


x 0 (t) = Ae at 


(5.4) 


From equations (5. 1) (5.4), it can be seen that the stability of a dynamic system 
depends upon the sign of the exponential index in the time response function, which 
is in fact a real root of the characteristic equation as explained in section 5.1.1. 
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Fig. 5.2 Graphical representation of stable and unstable time responses. 
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5.1.1 Stability and roots of the characteristic equation 

The characteristic equation was defined in section 3.6.2 for a second-order system as 

as 2 + bs + c — 0 (5.5) 



The roots of the characteristic equation given in equation (5.5) were shown in section 
3.6.2. to be 



s\, s 2 = 



— b ± Vb 2 — 4 ac 
2 a 



(5.6) 



These roots determine the transient response of the system and for a second-order 
system can be written as 

(a) Overdamping 



(b) Critical damping 



si = -a i 



S2 = -02 



(c) Underdamping 



= s 2 = —o 



(5.7) 



(5.8) 



si, s 2 = — <r± jw (5.9) 

If the coefficient b in equation (5.5) were to be negative, then the roots would be 

S\, s 2 = +a±)uj (5.10) 

The roots given in equation (5.9) provide a stable response of the form given in 
Figure 5.2(a) and equation (5.1), whereas the roots in equation (5.10) give an 
unstable response as represented by Figure 5.2(b) and equation (5.2). 

The only difference between the roots given in equation (5.9) and those in equation 
(5.10) is the sign of the real part. If the real part a is negative then the system is stable, 
but if it is positive, the system will be unstable. This holds true for systems of any 
order, so in general it can be stated: ‘If any of the roots of the characteristic equation 
have positive real parts, then the system will be unstable’. 



5.2 The Routh-Hurwitz stability criterion 

The work of Routh (1905) and Hurwitz (1875) gives a method of indicating the 
presence and number of unstable roots, but not their value. Consider the character- 
istic equation 



a„s" + a„-\f 1 + • • • + ci\s + oq = 0 



(5.11) 
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The Routh-Hurwitz stability criterion states: 

(a) For there to be no roots with positive real parts then there is a necessary, but not 
sufficient, condition that all coefficients in the characteristic equation have the 
same sign and that none are zero. 

If (a) above is satisfied, then the necessary and sufficient condition for stability is either 

(b) all the Hurwitz determinants of the polynomial are positive, or alternatively 

(c) all coefficients of the first column of Routh’s array have the same sign. The 
number of sign changes indicate the number of unstable roots. 

The Hurwitz determinants are 





D x 


= Cl\ 


d 2 = 


Cl\ 

00 


03 

02 














01 


03 


05 


07 


Cl\ 


03 


a 5 




00 


02 


04 


06 


00 


Cl2 


04 


d 4 = 




01 


03 


05 




a\ 


03 








02 


04 



(5.12) 



Routh’s array can be written in the form shown in Figure 5.3. 
In Routh’s array Figure 5.3 



l 


0«— 1 


0/2— 3 


Q-n— 1 


0« 


0/2— 2 


1 

Cl = — 


b\ 


t>2 


b\ 


0«— l 


0/2— 3 



bl = 



1 

Cl n —\ 



An— 1 An— 5 

Clfi 0^ — 4 



etc. 



c 2 = t~ 
b 1 



b\ b 3 

0/2— 1 0«— 5 



etc. 



(5.13) 

(5.14) 



Routh’s method is easy to apply and is usually used in preference to the Hurwitz 
technique. Note that the array can also be expressed in the reverse order, commen- 
cing with row s n . 



Pi 

<7i 



Ci 

Pi 

3 / 7-1 

a„ 



3 / 7-3 

a n -2 



3 / 7-5 

3 / 7-4 



Fig. 5.3 Routh's array. 
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Example 5.1 (See also Appendix 1, A1.5) 

Check the stability of the system which has the following characteristic equation 

s 4 + 2s 3 + s 2 + 4s + 2 = 0 (5.15) 

Test T. All coefficients are present and have the same sign. Proceed to Test 2, i.e. 
Routh’s array 



2 (5.16) 

4 

1 2 




The bottom two rows of the array in (5.16) are obtained from the characteristic 
equation. The remaining coefficients are given by 






2 4 
1 1 



= 2 (2 - 4) = -l 



bl= 2 



2 0 

1 2 



= -(4-0) = 2 



ci = -1 



b 3 = 0 



-1 2 
2 4 



= — 1(— 4 — 4) = 



c 2 = 0 



cl i = 



8 0 

-1 2 



= 4 ( 16 - 0 ) = 2 



(5.17) 

(5.18) 

(5.19) 

(5.20) 

(5.21) 

(5.22) 



In the array given in (5. 16) there are two sign changes in the column therefore there 
are two roots with positive real parts. Hence the system is unstable. 



5.2.1 Maximum value of the open-loop gain constant for the 
stability of a closed-loop system 



The closed-loop transfer function for a control system is given by equation (4.4) 



£ (s) = G{s) 

R ’ 1 + G(s)H(s) 



(5.23) 



In general, the characteristic equation is most easily formed by equating the denomi- 
nator of the transfer function to zero. Hence, from equation (5.23), the characteristic 
equation for a closed-loop control system is 

1 + G(s)H(s) = 0 



(5.24) 
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Proportional Control 

Controller Valve Plant 




Fig. 5.4 Closed-loop control system. 



Example 5.2 (See also Appendix 1, examp52.ni) 

Find the value of the proportional controller gain K\ to make the control system 
shown in Figure 5.4 just unstable. 

Solution 

The open-loop transfer function is 





(5.25) 


The open-loop gain constant is 




II 

oo 


(5.26) 


giving 




GWHW = ^/ i+ 2) 


(5.27) 


From equation (5.24) the characteristic equation is 




1 + , K = 0 
^(^ 2 + s + 2) 


(5.28) 


or 




,s(.r + ^ + 2) + A: = 0 


(5.29) 


which can be expressed as 




■s 3 + s 2 + 2s + K = 0 


(5.30) 



The characteristic equation can also be found from the closed-loop transfer function. 
Using equation (4.4) 

g(,) = 

R J I + (j(s)II(s) 



Given the open-loop transfer function in equation (5.27), where H(s) is unity, then 



C 

R 



CO 



K 

i(.5 2 +.S+2) 

1 I * 

' s(s 2 +.s+2) 



(5.31) 
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Multiplying numerator and denominator by s(s 2 + s + 2) 



C 

R 



(s) 



K 

5(5 2 + 5 + 2) + A' 



(5.32) 



-( S \ = 

s 3 + s 2 + 2s + K 

Equating the denominator of the closed-loop transfer function to zero 

5 3 + s 2 + 2s + K = 0 



(5.33) 



(5.34) 



Equations (5.30) and (5.34) are identical, and both are the characteristic equation. It 
will be noted that all terms are present and have the same sign (Routh’s first 
condition). Proceeding straight to Routh’s array 



where 



5° K 

s l (2 - A) 

s 2 1 K 

s 3 1 2 



b\ 



1 K 
1 2 



(2 - A) 



b 2 = 0 

Ci=K 



(5.35) 



To produce a sign change in the first column, 

K> 2 (5.36) 

Elence, from equation (5.26), to make the system just unstable 

K\ = 0.25 

Inserting (5.36) into (5.30) gives 

s 3 + s 2 + 2s + 2 = 0 

factorizing gives 

fr + 2 )(. s + 1) = 0 

hence the roots of the characteristic equation are 

5 = 0 ± jV2 

and the transient response is 

c(t) = Ae~‘ + B sin(\/2? + </>) (5.37) 

From equation (5.37) it can be seen that when the proportional controller gain K\ is 
set to 0.25, the system will oscillate continuously at a frequency of y/2 rad/s. 
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5.2.2 Special cases of the Routh array 

Case 1: A zero in the first column 

If there is a zero in the first column, then further calculation cannot normally proceed 
since it will involve dividing by zero. The problem is solved by replacing the zero with 
a small number e which can be assumed to be either positive or negative. When the 
array is complete, the signs of the elements in the first column are evaluated by 
allowing e to approach zero. 

Example 5.3 

s 4 + 2s 3 + 2s 2 + 4s + 3 = 0 

s° 3 

.v 1 4 — 6/e 

s 2 £ 3 (5.38) 

s 3 2 4 

s 4 1 2 3 

Irrespective of whether e is a small positive or negative number in array (5.38), there 
will be two sign changes in the first column. 



Case 2; All elements in a row are zero 

If all the elements of a particular row are zero, then they are replaced by the 
derivatives of an auxiliary polynomial, formed from the elements of the previous row. 



Example 5.4 

s 5 + 2s 4 + 6s 3 + 12s 2 + 8s + 16 = 0 

s° 16 

.v 1 8/3 

s 2 6 16 

s 3 8 24 

s 4 2 12 16 

s 5 16 8 



(5.39) 



The elements of the .y 3 row are zero in array (5.39). An auxiliary polynomial P(s) is 

therefore formed from the elements of the previous row (/). 

i.e. 

Pis) = 2 s 4 + 12s 2 + 16 

= 8 s 3 + 24.y (5.40) 

The coefficients of equation (5.40) become the elements of the ,s 3 row, allowing the 
array to be completed. 
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5.3 Root-locus analysis 
5.3.1 System poles and zeros 

The closed-loop transfer function for any feedback control system may be written in 
the factored form given in equation (5.41) 

C G'(.v) _ K c js - z c \ )(.? - z c2 ) ...(s- z c „ ) 

R ~ 1 + G(s)H(s) ~ {s - p c i){s - p c2 ) ...is- Pen) ^ ’ 

where 5 = p c \, pa, p cn are closed-loop poles, so called since their values make 
equation (5.41) infinite (Note that they are also the roots of the characteristic 
equation) and s = z c i, z c2 , . . . , z C n are closed-loop zeros, since their values make 
equation (5.41) zero. 

The position of the closed-loop poles in the .v-plane determine the nature of the 
transient behaviour of the system as can be seen in Figure 5.5. Also, the open-loop 
transfer function may be expressed as 

= (5.42) 

is - poi)is - p 02 ) . . . is - p Qn ) 

where z 0 i, zq 2 , ■ ■ ■ , -o « are open-loop zeros and p 0 1 , p 02 , . . . , po„ are open-loop poles. 




Fig. 5.5 Effect of closed-loop pole position in the s-plane on system transient response. 
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5.3.2 The root locus method 

This is a control system design technique developed by W.R. Evans (1948) that 
determines the roots of the characteristic equation (closed-loop poles) when the 
open-loop gain-constant K is increased from zero to infinity. 

The locus of the roots, or closed-loop poles are plotted in the s-plane. This is a 
complex plane, since s = a ± jw. It is important to remember that the real part a is 
the index in the exponential term of the time response, and if positive will make the 
system unstable. Hence, any locus in the right-hand side of the plane represents an 
unstable system. The imaginary part co is the frequency of transient oscillation. 

When a locus crosses the imaginary axis, a = 0. This is the condition of marginal 
stability, i.e. the control system is on the verge of instability, where transient oscilla- 
tions neither increase, nor decay, but remain at a constant value. 

The design method requires the closed-loop poles to be plotted in the s-plane as K 
is varied from zero to infinity, and then a value of K selected to provide the necessary 
transient response as required by the performance specification. The loci always 
commence at open-loop poles (denoted by x) and terminate at open-loop zeros 
(denoted by o) when they exist. 

Example 5.5 

Construct the root-locus diagram for the first-order control system shown in 
Figure 5.6. 

Solution 

Open-loop transfer function 

G(s)H{s) = K (5.43) 

Ts 

Open-loop poles 

s = 0 

Open-loop zeros: none 

Characteristic equation 

1 + G(s)H(s ) = 0 

Substituting equation (5.3) gives 

<4=o 

i.e. Ts+K=0 (5.44) 




Fig. 5.6 First-order control system. 
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Fig. 5.7 Root-locus diagram for a first-order system. 



Roots of characteristic equation 




(5.45) 



When K is varied from zero to infinity the locus commences at the open-loop pole 
j = 0 and terminates at minus infinity on the real axis as shown in Figure 5.7. 

From Figure 5.7 it can be seen that the system becomes more responsive as K is 
increased. In practice, there is an upper limit for K as signals and control elements 
saturate. 



Example 5.6 

Construct the root-locus diagram for the second-order control system shown in 
Figure 5.8. 

Open-loop transfer function 

G(s)H(s) = — (5.46) 
s(s + 4) 



Open-loop poles 



Open-loop zeros: none 

Characteristic equation 



s = 0,-4 



1 + G(s)H(s) = 0 




Fig. 5.8 Second-order control system. 
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Table 5.1 Roots of second-order characteristic 
equation for different values of K 



K 


Characteristic equation 


Roots 


0 


s 2 + 4.s = 0 


^T 

1 

o' 

II 

03 


4 


s 2 + 4s + 4 = 0 


s= —2 ± jO 


8 


j 2 + 4s + 8 = 0 


.? = — 2±j2 


16 


s 2 + 4.s+ 16 = 0 


s = -2±j3.46 



Substituting equation (5.4) gives 



i.e. 



1 + 



K 



s(s + 4) 



= 0 



4s + K = 0 



(5.47) 



Table 5.1 shows how equation (5.7) can be used to calculate the roots of the 
characteristic equation for different values of K. Figure 5.9 shows the corresponding 
root-locus diagram. 

In Figure 5.9, note that the loci commences at the open-loop poles ( s = 0, —4) 
when K = 0. At K = 4 they branch into the complex space. This is called a break- 
away point and corresponds to critical damping. 




Fig. 5.9 Root locus diagram for a second-order system. 



